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I.  Introduction 


To  produce  a  high  performance  ultrasonic  transducer  requires  the  transducer  material,  which 
performs  the  energy  conversion  between  the  mechanical  form  and  electrical  form,  to  have  a  high 
electromechanical  coupling  factor,  broad  operation  frequency  bandwidth,  and  adjustable  acoustic 
impedance  which  can  be  tuned  to  match  that  of  the  medium.  With  single  phase  piezoelectric  materials, 
it  is  difficult  to  simultaneously  meet  all  the  requirements.  Piezoceramic  polymer  composite  materials, 
which  combine  the  high  electromechanical  activity  of  piezoceramics  and  the  low  acoustic  impedance  of 
polymeric  materials,  have  provided  new  opportunities  to  meet  these  requirements.1,2  Since  their 
inception  in  the  seventies,  the  piezoceramic  polymer  composites  have  become  one  of  the  most 
important  transducer  materials  and  are  being  widely  used  in  many  areas  such  as  medical  imaging, 
mondestructive  evaluation  of  materials,  and  underwater  vision,  etc. 


Being  a  diphasic  material,  the  properties  of  a  piezo-composite  can  be  tailored  over  a  wide  range 
by  adjusting  the  material  properties  and  geometric  shapes  of  constituent  phases.3’4  It  has  also  been 
observed  that  the  properties  of  a  composite  vary  with  frequency.5  The  challenge  of  understanding  the 
seemingly  complex  relationship  between  the  performance  of  a  composite  and  the  properties  of  its 
constituents  and  the  great  opportunities  provided  by  these  materials  have  stimulated,  in  the  past  two 
decades,  extensive  investigations,  both  experimental  and  theoretical,  on  this  class  of  materials. 

The  classic  work  of  Newnham  et  al.,6  which  classified  piezocomposite  materials  according  to 
the  connectivity  of  the  constituent  phases,  has  greatly  facilitated  the  analysis  of  the  composites  as  the 
connectivity  is  one  of  the  key  parameters  in  determining  the  performance  of  a  composite.  For 
composites  with  1-3  and  2-2  connectivities,  both  analytical  and  finite  element  modeling  have  been 
earned  out  which  have  provided  useful  guidelines  in  the  design  of  composite  transducers.7'13  The 
isostrain  models  developed  by  Smith  et  al.3  and  Hashimoto  et  al.4  linked  the  material  parameters  of  the 
constituents  to  the  effective  piezoelectric  properties  of  1-3  and  2-2  composites  respectively  and 
predicted  that  the  thickness  coupling  factor  lq  of  a  composite  can  approach  the  longitudinal  coupling 
factor  £33  of  the  piezoceramic  rod  and  £33  of  the  plate  respectively,  which  is  in  good  agreement  with 
experiments  for  composites  with  a  high  aspect  ratio  t/d,  where  t  is  the  thickness  and  d  is  the 
periodicity  of  the  composites.  Auld  et  al.,  using  the  Roquet  theory,  investigated  wave  propagation  in 
both  the  2-2  and  1-3  composites  and  showed  that  due  to  the  periodic  structure  of  these  composites, 
there  exist  pass  bands  and  stop  bands,  similar  to  the  band  structure  in  a  crystal  solid,  and  that  there  are 
piezoelectric  resonances  associated  with  the  stop  band  edge  resonances.710  For  the  design  of 
composite  transducers,  the  recognition  of  the  existence  of  these  modes  and  the  precise  prediction  of 
their  frequencies  are  of  prime  importance  since  quite  often,  it  is  the  interference  of  these  modes  with 
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the  thickness  mode  of  a  composite  that  deteriorates  the  performance,  especially  at  high  frequency 
operations.  Craciun  et  al.  examined  the  coupling  between  these  lateral  modes  and  thickness  mode  using 
a  phenomenological  approach  and  the  results  provided  qualitative  understanding  between  the  coupling 
of  the  two  modes  and  the  material  properties.11  The  results  from  these  investigations  have  played 
important  roles  in  the  development  of  ultrasonic  composite  transducers.  However,  due  to  the 
approximations  used  in  the  analysis,  there  are  severe  limitations.  For  instance,  various  features  related 
to  the  dynamic  behavior  of  a  composite  transducer  were  not  treated  in  a  consistent  manner,  and  the 
effect  of  finite  thickness  of  a  composite  on  the  material  properties,  that  has  been  shown  to  be  crucial  in 
determining  the  performance  of  a  piezocomposite,  cannot  be  treated  in  these  analyses.  To  address  these 
realistic  issues  of  a  composite  material,  finite  element  analysis  has  been  employed  by  many  authors.12" 
15  For  example,  the  dispersion  curves  have  been  evaluated  for  various  modes  in  a  composite  and  the 
^dependence  of  the  electromechanical  coupling  factor  on  the  ceramic  volume  content  and  the  ceramic  rod 
shape  was  investigated. 

For  a  piezoceramic  polymer  composite,  it  has  to  be  recognized  that  it  is  the  ceramic  phase 
which  performs  the  energy  conversion  between  the  electric  and  mechanical  forms  and  the  polymer 
phase  merely  acts  as  a  carrier  which  transfers  acoustic  energy  between  the  piezoceramic  and  the 
external  medium.  Hence,  if  the  elastic  coupling  between  the  two  constituents  is  not  very  effective,  even 
if  the  material  exhibits  a  perfect  acoustic  impedance  matching  with  the  medium  and  the 
electromechanical  coupling  factor  is  large,  the  electromechanical  performance  of  the  material  is  still 
poor.  These  observations  clearly  indicate  that  in  modeling  piezoceramic  polymer  composites,  one 
cannot  simply  use  an  effective  medium  approach  and  has  to  take  into  account  explicitly  this  internal 
degree  of  freedom. 

The  objective  of  this  program  is  to  develop  a  realistic  model  which  can  address  these  practical 
issues  and  to  conduct  systematic  experimental  study  to  verify  the  model  results  and  to  provide  new 
insight  into  composite  materials  for  both  underwater  and  ultrasonic  transducer  applications. 

Through  this  program,  an  analytical  model  on  the  dynamic  problem  of  a  piezocomposite 
material  with  the  2-2  connectivity  and  1-3  connectivity  has  been  developed.16'18  In  the  model,  we 
avoided  the  approximations  made  in  the  earlier  works  and  hence,  can  address  the  dynamic  responses 
of  2-2  piezocomposites,  such  as  the  frequencies  of  various  modes,  the  mode  coupling,  the 
electromechanical  coupling  factor,  the  vibration  profiles  of  composites  under  different  external  driving 
conditions,  etc.,  in  a  realistic  and  consistent  manner. 
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Figure  1.  (a)  Schematic  drawing  of  a  2-2  piezoceramic  polymer  composite  which  is  unbounded  in 
aU  three  orthogonal  directions.  The  period  of  the  composite  is  d.  (b)  Schematic  drawing  of  a  2-2 
piezoceramic  polymer  composite  plate  with  a  thickness  t  situated  in  air.  An  external  voltage  is 
applied  to  the  composite.  ~ 


Figure  2.  The  first  and  second  branches  of  the  dispersion  curves  (solid  lines)  and  the  experimentally 
measured  thickness  mode  (solid  circles)  and  lateral  mode  (open  circles)  for  a  2-2  composite  made  of 
PZT-5H  and  Spurr  epoxy  with  44%  ceramic  content.  After  the  crossover  region,  the  thickness  mode 

to  the  second  branch.  The  criterion  for  the  thickness  mode  is  that  it  is  the  one  with  the  larger 
effective  coupling  factor. 
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For  example,  one  of  the  misconceptions  in  the  early  studies  of  the  dynamic  behavior  of 
piezocomposites  is  the  direct  linkage  between  the  non-uniform  surface  vibration  profile  in  a  composite 
and  the  aspect  ratio  of  the  ceramic  plate  (in  2-2  composites)  or  rod  (in  1-3  composites).  Based  on  the 
model  results,  it  will  be  shown,  which  was  verified  by  experiment,  that  for  a  composite  plate,  as  long 
as  the  thickness  resonance  frequency  is  below  that  of  the  lateral  mode,  there  is  always  a  frequency  near 
the  thickness  resonance  where  the  vibration  profile  of  the  composite  is  uniform.  The  influence  of  the 
aspect  ratio  is  on  the  frequency  bandwidth  in  which  the  polymer  and  ceramic  vibrate  in  unison. 

Our  results  also  show  that  there  exist  a  series  of  modes  associated  with  the  periodic  structure  of 
the  composite,  which  is  beyond  the  stop-band  edge  resonance  prediction.  In  1-3  composites,  the 
experimental  results  indicate  that  the  nature  of  these  modes  is  far  more  complicated  than  the  simple 
‘picture  of  stop-band  edge  resonance.  Further  theoretical  modeling  is  clearly  needed  in  this  area  since 
these  modes  play  a  crucial  role  for  a  broadband  transducer. 

The  approach  taken  here  to  solve  the  vibration  problem  in  a  finite  thickness  composite  plate  is 
based  on  the  method  of  partial  wave  expansion  where  various  elastic  and  electric  fields  are  expended  in 
terms  of  the  eigenmodes  of  the  structure  as  shown  in  figure  1(a)  and  the  coefficients  for  each  eigen¬ 
mode  are  determined  by  the  boundary  conditions  at  x3  =  t  /2  (figure  1(b)). 

In  this  program,  we  also  examined  how  the  material  properties  at  the  two  sides  of  the 
interface,  i.e.,  the  properties  of  a  piezoceramic  polymer  composite  and  a  uniform  medium,  either  a 
liquid  or  a  solid,  influence  the  wave  transmission  and  reflection  at  the  interface  and  the  acoustic  energy 
transfer  inside  the  composite.  These  are  directly  related  to  the  design  of  the  matching  layer  and 
backing  for  a  composite  transducer.  From  the  results,  it  is  shown  that  the  thickness  of  the  matching 
layer  should  be  less  than  the  conventional  X/4.  It  is  also  shown  that  the  input  acoustic  impedance  of  a 
composite  depends  strongly  on  the  frequency  and  the  elastic  properties  of  the  medium  which  the 
composite  is  in  contact  with. 

This  report  is  organized  as  the  following.  In  the  section  n,  the  results  on  the  dynamic  behavior 
of  a  2-2  composite  transducer  are  presented.  We  will  first  present  the  details  of  derivation  of  the 
eigenfunctions  for  a  2-2  piezoceramic  polymer  composite.  Based  on  the  results,  the  wave  propagation 
in  an  unbounded  2-2  composite  is  analyzed.  For  finite  thickness  composites,  various  features  of  a 
piezocomposite  plate  under  an  external  driving  electric  field  in  both  air  and  water  media  are  treated. 
The  response  behavior  of  the  composite  plate  under  external  pressure  is  also  treated.  Experiments 
were  conducted  to  provide  comparison  with  the  theoretical  results.  From  the  results,  we  will  show 
explicitly  how  the  aspect  ratio  t/d  influences  the  transduction  performance  of  a  composite  plate. 
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In  the  section  m,  the  acoustic  impedance  and  the  wave  propagation  at  the  interface  between  a 
2-2  composite  and  an  acoustic  medium,  either  a  solid  or  a  fluid  are  analyzed.  The  general  formula  used 
to  treat  the  wave  propagation  at  the  interface  is  briefly  presented  and  the  results  for  the  fluid-composite 
interface  are  derived.  The  definition  of  the  input  acoustic  impedance  of  the  composite  at  the  interface  is 
also  discussed.  The  wave  propagation  at  a  interface  between  a  solid  and  a  piezocomposite  is  analyzed. 
The  issues  related  to  the  anti-reflection  quarter  wave  matching  layer  at  a  composite-medium  interface 
are  examined.  And  a  brief  summary  and  discussion  are  presented  at  the  end  of  the  section. 

In  the  section  IV,  the  results  on  1-3  composites  are  presented.  The  effect  of  the  loss  in  the 
elastic  properties  in  the  polymer  phase  is  considered  and  it  was  found  that  due  to  the  acoustic  loss,  the 
^dispersion  curve  for  the  lateral  mode  in  the  composite  shows  marked  change.  The  experimental  results 
related  to  the  lateral  modes  are  also  presented. 

A  summary  will  be  presented  on  the  section  V. 


n.  Dynamic  Behavior  of  a  Finite  Thickness  2-2  Piezoceramic -Polymer  Composite  Plate: 
2T — General  Solution  of  the  Waves  in  a  2-2  Piezocomposite: 


Shown  in  figure  1(a)  is  a  schematic  drawing  of  a  2-2  composite  (unbounded  in  the  x3- 
direction),  where  plates  of  piezoceramic  and  polymer  form  a  parallel  array.  The  coordinate  system  is 
chosen  such  that  the  x3-axis  is  along  the  ceramic  poling  direction,  the  x,-axis  is  perpendicular  to  the 
ceramic  polymer  interface,  and  x2-axis  is  in  the  plane  of  the  plates.  For  a  typical  2-2  composite,  the 
dimensions  in  the  x,-  and  x2-directions  are  much  larger  than  the  period  d  and  thickness  t.  In  the 
treatment  here,  they  can  be  taken  as  infinite  without  much  error  in  the  results.  Under  these  conditions, 
the  composite  is  clamped  in  the  x2-direction  so  that  S2=0,  where  S2  is  the  strain  in  the  x2-direction,  and 
the  problem  becomes  a  two  dimensional  one  with  no  dependence  on  the  x2-coordinate. 

The  governing  equations  for  the  dynamics  of  a  2-2  composite  are19'21 
|  fj%  d2“i 
&1  <&3  P  dt2 

(1) 

<&!  <&3  P  dt 2 
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The  symbols  adopted  in  this  paper  are  summarized  here:  T,  and  S;  are  the  stress  and  strain  tensor 
components,  where  the  Voigt  notation  is  used,  u;  is  the  elastic  displacement  vector,  p  is  the  density, 
and  D;  is  the  electric  displacement  vector,  E;  is  the  electric  field.  The  relevant  material  coefficients  are: 
eij  is  the  piezoelectric  coefficient,  cy  is  the  elastic  stiffness,  and  £{  is  the  dielectric  permittivity.  Equation 
(1)  holds  for  both  the  polymer  and  ceramic  phases. 

The  constitutive  equations,  relating  stress  T,  strain  S,  electric  displacement  D,  and  electric  field 

E,  are 

m  =  Ws]-kl£]  (2a) 

_  [D)  =  [e][S]+[es][£]  (2b) 

For  the  polymer  phase,  e^  in  the  equation  (2)  are  zero.  The  superscripts  E  and  S  indicate  that  the 
coefficients  are  under  the  constant  E  field  and  constant  strain  conditions,  respectively.  Here  [et]  is  the 

transposed  [e]  array.  Under  the  quasi-electrostatic  approximation,  the  electric  field  E  can  be  expressed 
as 

E  =  -VO  ( 3 ) 

where  <E>  is  the  electrical  potential. 

Combining  eqs.  (1),  (2),  and  (3)  yields  three  second  order  differential  equations,  governing  the 
elastic  displacement  Up  u3,  and  electric  potential  <I>  in  the  ceramic  plate,  respectively: 

4M3,1 1  +  (4  +  4  )«U3  +  4«3.33  +  (*33^.33  +  *15®.  1 1 )  =  P«3 

CUUU 1  +  (4  +  4  )«3.13  +  4«U33  +  («31  +  el5  )^,13  =  P«i  (  4  ) 

ei5U3,U  +  (*15  +  e3lKl3  +  e33U3  33  -  (£33^.33  +  4^.11)  =  0 

For  the  polymer  phase,  e;j-  in  eq.  (4)  should  be  taken  to  zero. 

For  an  unbounded  composite,  the  solutions  to  eq.  (4)  have  the  form 
«3  =  A  &xp(j(hxl  +  (3x3  -  cot)) 

ui  =  BQxp(j(hxl  +  px3  -  cot))  ( 5  ) 

O  =  C  exp(j(hxl  +  px3  -  cot)) 

where  A,  B,  and  C  are  three  constants,  co  is  the  angular  frequency,  h  and  (3  are  the  wave  vector 
components  in  the  x{-  and  x3-  directions,  respectively. 


Substituting  eq.  (5)  into  (4)  yields  three  homogeneous  equations  with  the  undetermined 
constants  A,  B,  and  C, 

MM  =  0  (  6  ) 


(7) 


where  [AJ=[A,  B,  C]T  and 

(4P2+cEuh2-pco2  (cft  +  c&hp  e3ip2+el5h2  ' 

[Mij]=  (4+4)^  cEh2  +  c^p2-pco2  (el5+e3l)h{3 

l  exP2+ei5h2  (eis+en  )hP  -(ef,  h2+£s33P2)j 

The  condition  for  a  nontrivial  solution  is  such  that  the  determinant  of  the  coefficient  matrix  vanishes, 
i.e., 

m=o  (8) 

Equation  (8)  is  a  cubic  equation  of  h2.  For  a  given  co  and  (3 ,  eq.  (8),  in  general,  has  three  roots  of  h2, 
denoted  as  hj,  h3  and  h3 ,  corresponding  to  the  quasi-electromagnetic,  quasi-longitudinal,  and  quasi¬ 
shear  waves  in  the  piezoelectric  plate,  respectively.  For  each  h',  the  ratio  among  A,  B,  and  C  can  be 
determined  from  eq.  (6).  Since  we  are  concerned  only  with  the  waves  in  the  x3  direction  which 
Correspond  to  piezo-active  modes  in  a  finite  thickness  composite  plate,  the  general  solutions  become 
*fj  =  X  R‘  fi  C0s<Ac*i )  sin(/?x3 ) 

i 

4  =  X  sin0fo)cos(/fc3 )  ( 9 ) 

f 

=  X^  cos(h'xx)sm(px3) 

i 

where  i  runs  from  1  to  3.  fi(  g,  and  t;  are  the  cofactors  of  Akl(i),  Ak2(i),  and  Au(i)  of  the  determinant 

(8)  (where  h  is  replaced  by  h°  for  i=l,  2,  and  3,  respectively).  The  exp(-jcot)  term  in  eq.  (9)  is 
omitted. 

Following  the  similar  procedure,  the  solutions  for  the  polymer  phase  can  be  obtained  (the 
center  of  the  polymer  plate  is  at  x,  =  d/2) 

«f  =  IX fiPcos(hE(xl  -  j))sin (flx3) 

4  =  X  PPSi  sin(4Ui  ~  —))  cos (fix3 )  ( 10 ) 

i  ^ 

Op  =  Cp cosh(/J(x,  -^))sin(^x3) 

where  i  =1,  2.  fp  and  gp  are  the  cofactors  of  Akl(i),  A^fi)  of  the  determinant  (8)  with  all  the  material 
parameters  replaced  by  those  of  the  polymer  phase,  and 

(h1}1)2  =  (k£)2  -  p2  and  (h|)2  =  (k£)2  -  p2 

where  kL  =  — ,  kT  =  — ,  vpL  and  vp  are  the  longitudinal  and  shear  wave  velocities  of  the  polymer 

VL  Vr 

phase,  respectively. 
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The  expressions  of  the  stresses  and  the  electric  displacement  in  the  ceramic  plate  can  be 
obtained  by  substituting  eq.  (9)  into  eq.  (2): 

T\  ~  X Tn  (O^f  cos(h‘:.xl)cos(jix3 ) 

i 

T3C  =  X  ^3  (Otf  COS (hfxl  )COS(pXi)  (  1 1  ) 

i 

if  =  sinfAfxJsinf/ir.) 

f 


where 

(0  =  cuhfgf  +  +  e3  xtf)p 

where  i=l,  2,  3  for  ceramic  phase.  The  electric  displacements  are: 
D\  =  (i)l^  sinQif  x{)sm(px3) 

i 

D3  =  X  AC(0  cos(hfx:)cos({3x3) 


(12) 


(13a) 

(13b) 


where  Dp  (i)  =  -el50gP  +  hf, ftc)  +  e5nh?tf  andZJ,c(i)  =  e3lhfgf  +  -  e&  ^.Similar  expressions 

can  be  obtained  for  the  polymer  plate.  For  instance,  the  electric  displacements  are: 

Di  =  ~euPCP sinh(/?(jCj  -  j))sin (J3x3)  (  14C ) 

D3  =  -*\\PC?  cosh(j8(x,  -  j))  cos(fix3 )  ( 14d ) 


Making  use  of  the  boundary  conditions  at  the  ceramic  polymer  interface  (x(=  vd/2,  where  v  is 
the  volume  fraction  of  the  ceramic  in  a  composite),  which  are 

=  7;C  =  7;P  and  7/  =  7/  (15a) 

and  Oc  =  <DP  and  D,c  =  £>f  ( 15b  ^ 

Rj  and  Cp  (i  =1,  2,  3  and  j=l,  2)  can  be  determined.  Expanding  Eqs.  (15)  in  terms  of  eqs.  (9), 
(10),  (11),  (13)  and  (14)  yields  the  following  six  homogeneous  linear  equations: 
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(16) 


f*i, 

*12 

*13 

*14 

*15 

0  > 

(R?) 

*2. 

■*22 

*23 

*24 

*25 

0 

*3. 

*32 

*33 

*34 

*35 

0 

*,c 

*41 

*42 

*43 

*44 

*45 

0 

rz 

0 

0 

*53 

*54 

*55 

*56 

*3C 

,0 

0 

*63 

*64 

*65 

*66  ) 

c) 

where  the  matrix  elements  Ky  are  functions  of  (3  ,a>,  d,  v,  and  the  material  parameters  of  both  the 
polymer  and  piezoceramic.  The  condition  that  the  determinant 

I  Ky  |  =0  (17) 

yields  the  relationship  between  0)  and  (3,  the  dispersion  relations  in  the  composite.  For  each  pair  of  © 
^_and  (3  ,  the  relationships  among  Rf ,  R?  and  Cp  can  be  obtained.  Hence,  the  various  stress,  strain, 
electric  field  distributions  in  the  composite  can  be  determined. 

~~~~  Dispersion  Curves.  Modes  and  Mode  Coupling  of  a  2-2  Piezocomposite! 

Eq.  (17)  allows  us  to  determine  the  relationship  between  (3  and  ©,  the  dispersion  curves,  for  a 
composite  if  the  materials  parameters  of  piezoceramic  and  polymer  and  the  geometric  parameters,  such 
as  d  and  v,  are  known.  Eq.  (17)  is  a  transcendental  function  which  cannot  be  solved  analytically.  A 
computer  code  was  developed  and  the  dispersion  curves  were  evaluated  numerically. 

Shown  in  figure  2  (the  solid  curves)  are  the  two  lowest  branches  of  the  dispersion  curves  for  a 
2-2  composite  made  of  PZT-5H  piezoceramic  and  Spurr  epoxy  with  the  ceramic  volume  fraction  of 
44%.  For  all  the  composites  discussed  in  this  paper,  except  otherwise  specified,  PZT-5H 
piezoceramic  and  Spurr  epoxy  are  used  as  the  constituents.  The  parameters  of  PZT-5H  and  Spun- 
epoxy  are  presented  in  Table  I.  The  general  trend  of  the  dispersion  curves  in  figure  2  resembles  that  of 
the  symmetric  Lamb  waves  in  a  plate.19  As  will  be  shown  later,  at  small  Pd  limit,  the  first  branch 
corresponds  to  the  longitudinal  wave  propagation  along  the  x3-direction,  that  is,  u3  is  the  dominant 
displacement  which  is  more  or  less  uniform  and  in  phase  in  the  x, -direction,  and  its  phase  velocity  is 
the  effective  longitudinal  wave  velocity  of  the  composite.  The  second  branch  corresponds  to  the  lateral 
resonance  which  arises  from  the  periodic  structure  of  the  composite,  that  is,  it  is  a  stationary  shear 
wave  along  the  x,-direction,  and  it  can  be  shown  that  the  displacement  u3  of  the  polymer  phase  is  much 
larger  than  that  of  the  ceramic  phase  and  the  phase  difference  between  them  is  180°.  Hence,  it  can  be 
related  to  the  so  called  stop-band  edge  resonance  as  predicted  by  Auld  et  al.7-8 
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Table  I.  The  material  properties  of  PZT-5H  and  Spurr  epoxy  used  for  the  2-2  composites  in  the 
investigation. 

PZT-5H:  e33  =  23.09  C/m2,  e31  =  -  6.603  C/m2,  eI5  =  17.0  C/m2,  c„  =  12.72*  1010  N/m2,  = 

2.3*  1010  N/m2,  c33  =  11.74*10’°  N/m2,  cn  =  8.47*1010  N/m2,  Kn  =  1700,  K33  =  1470,  p  =  7500 
kg/m3. 

Spurr  epoxy:  cn  =  5.4*109  N/m2,  =  1.3*  109  N/m2,  p  =  1100  kg/m3. 


Under  the  assumption  that  the  wave  length  X  in  the  x3-direction  ((3  =2jt/  X)  is  equal  to  two 
times  the  composite  thickness  t  (the  condition  for  the  thickness  resonance),  the  theoretical  dispersion 
curves  can  be  compared  with  the  experimental  results,  obtained  from  composites  with  different 
thickness.  At  low  (3d  on  the  dispersion  curves,  where  d  is  the  period  of  the  composite,  the  frequency 
of  the  fundamental  thickness  mode  from  the  experiment  falls  on  the  first  branch  as  marked  by  the  black 
‘dots,  and  the  frequency  of  the  first  lateral  mode  falls  on  the  second  branch  as  marked  by  the  open 
circles.  After  the  crossover  region  B,  the  modes  interchange  the  positions  on  the  dispersion  curves 
where  the  thickness  mode  is  in  the  second  branch  while  there  is  a  weak  resonance  at  a  frequency  below 
the  thickness  mode  and  it  falls  on  the  first  branch.  Here,  the  thickness  mode  is  defined  as  the  one  with 
higher  electromechanical  coupling  factor.  The  mode  on  the  first  branch  gradually  diminishes  at  high  (3d 
values.  The  result  is  presented  in  figure  2.  Clearly,  there  is  an  excellent  accord  between  the 
theoretically  derived  resonance  frequencies  and  experimentally  observed  ones. 

Shown  in  figures  3(a)  and  figure  3(b)  are  the  3-dimensional  dispersion  curves  for  2-2 
composites  with  15%  and  44%  ceramic  volume  fraction,  respectively.  Obviously,  there  are  imaginary 
and  complex  branches  of  the  dispersion  curves.  The  modes  on  these  branches  are  non-propagating 
modes  (imaginary  branches)  and  attenuated  modes  (complex  branches),  respectively.  These  modes  do 
not  exist  in  an  unbounded  composite,  however,  they  are  important  in  the  vibration  problems  of  finite 
thickness  plates  as  well  as  in  semi-infinite  mediums,  where  they  correspond  to  the  evanescent  waves  at 
the  surface  which  will  be  discussed  later  in  the  paper. 

Figure  4  presents  the  effect  of  the  stiffness  of  the  polymer  phase  on  the  dispersion  curves  of  a 
composite  with  30%  ceramic  content  where  PZT-5H  is  used  as  the  piezoceramic  phase.  It  is 
apparent  that  the  frequency  position  of  the  second  branch  is  very  sensitive  to  properties  of  the  polymer 
phase.  This  is  quite  understandable  since  as  will  be  shown  later,  the  frequency  position  of  this  branch 
is  directly  related  to  the  shear  velocity  of  the  polymer  matrix  for  composites  at  this  ceramic  content. 
The  lower  the  shear  wave  velocity  of  the  polymer  phase  is,  the  lower  the  frequency  of  the  second 
branch  will  be.  Therefore,  for  a  transducer  operated  at  high  frequencies,  in  order  to  avoid  the 
interference  from  the  lateral  modes,  a  polymer  matrix  with  a  high  shear  velocity  should  be  utilized  even 
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Figure  4.  The  effect  of  the  properties  of  the  polymer  matrix  on  the  dispersion  curves  of  2-2 
composites  with  PZT-5H  ceramic,  derived  from  the  model.  The  properties  of  the  polymer  phase  are: 
curve  .  cn-  7.72*10  N/nr,  c44=  1.588*109  N/m2,  p  =  1160  kg/m3;  curve  2:  c,.=  1.3634*1010 
i /.m2’  C44=  3-432*10  N/m2,  p  =  1610  kg/m3;  curve  3:  cu=  3.173*109  N/m2,  cu-  0.696*  109 
N/m  ,  p  =  1060  kg/m  ;  curve  4:  cn=  1.622*109  N/m2,  C44=  1.646*10®  N/m2,  p  =  890  kg/m3. 
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though  the  thickness  coupling  factor  k,  of  this  composite  may  be  reduced  as  a  result  of  the  stiffer 
polymer  matrix.3,4 

In  order  to  further  elucidate  the  origins  of  these  resonant  modes,  it  is  instructive  to  make  a 
comparison  between  the  dispersion  curves  of  a  composite  and  those  of  single  piezoceramic  and 
polymer  plates  with  appropriate  boundary  conditions.23  In  general,  in  piezoceramic  polymer 
composites,  the  elastic  stiffness  of  the  piezoceramic  is  more  than  one  order  of  magnitude  higher  than 
that  of  the  polymer.  It  is  reasonable  to  assume  that  the  piezoceramic  plates  in  a  2-2  composite  are  stress 
free  at  the  ceramic-polymer  interface,  that  is,  the  stress  components  T,  and  T5  are  zero  at  the  interface. 
Similarly,  the  polymer  plates  can  be  approximated  as  under  the  fixed  boundary  conditions,  that  is,  u, 
and  u3  are  zero  at  the  interface.  Under  these  assumptions,  the  dispersion  curves  of  both  single  ceramic 
jmd  single  polymer  plates  are  calculated  for  the  plates  with  different  widths  corresponding  to  2-2 
composites  with  different  ceramic  volume  fractions.  The  results  are  presented  in  figures  5(a),  5(b),  and 
5(c)  corresponding  to  2-2  composites  with  15%,  44%  and  80%  ceramic  content,  respectively.  In  these 
figures,  the  solid  lines  are  the  dispersion  curves  of  the  2-2  composites,  the  dashed  lines  are  those  of 
the  ceramic  plate,  and  doted  lines  are  the  dispersion  curves  of  the  polymer  plate. 


The  results  from  these  figures  reveal  that  there  are  many  resemblance  between  the  dispersion 
curves  of  2-2  composites  and  the  dispersion  curves  of  the  single  ceramic  and  single  polymer  plates 
with  appropriate  boundary  conditions.  For  instance,  at  small  (3d,  the  first  branch  in  the  dispersion 
curves  of  the  2-2  composites  with  44%  and  80%  ceramic  volume  content  is  very  close  to  the  first 
branch  of  the  ceramic  plate.  On  the  other  hand,  for  2-2  composites  with  low  ceramic  volume  content,  a 
large  difference  between  these  two  is  found  for  the  first  branch.  These  results  are  consistent  with 
those  from  Smith  et  al.  and  from  Hashimoto  et  al.  based  on  an  effective  medium  model.3,4  For  the 
second  branch  which  corresponds  to  the  lateral  resonant  mode  in  a  composite  plate,  at  low  and  medium 
ceramic  volume  content,  it  is  close  to  the  first  branch  of  the  polymer  plate  which  frequency  at  small  pd 
is  equal  to  VT  /2 dp  where  V-f  and  dp  are  shear  wave  velocity  and  width  of  the  polymer  plate, 
respectively.  However,  for  composites  with  high  ceramic  volume  content  such  as  the  one  shown  here 
(80%),  the  second  branch  is  related  to  the  longitudinal  resonance  of  the  ceramic  plate  along  the  width 
(or  xr)  direction  while  the  shear  resonance  associated  with  the  polymer  plate  lies  on  the  third  branch  of 
the  dispersion  curves. 

Obviously,  the  coupling  between  the  two  phases  through  the  interface  boundary  conditions  will 
influence  the  dispersion  curves  of  the  waves  in  the  two  phases.  It  is  well  known  that  the  dispersion 
curves  for  the  uncoupled  waves  are  split  at  their  crossover  points  when  coupling  is  introduced.11  Far 
from  the  crossover  region,  the  coupled  wave  dispersion  curves  should  nearly  coincide  with  those  of 
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the  uncoupled  waves.  If  the  coupling  is  very  strong  as  for  the  15%  and  44%  piezocomposites,  the 
coupled  waves  exhibit  large  departure  from  the  uncoupled  curves  in  the  crossover  region,  which  is 
clearly  shown  in  figure  5(a)  and  figure  5(b).  While  for  the  composite  with  80%  ceramic  content,  the 
coupling  between  the  first  and  the  second  branches  is  through  the  coupling  of  P  wave  and  SV  wave  in 
the  ceramic  plate  where  the  interface  does  not  have  a  significant  effect  on  it. 

2.3.  Vibration  of  a  Finite  Thickness  Composite  Plate  under  an  F.lerthc  Field  in  Air 

In  the  previous  section,  the  properties  of  guided  wave  propagation  in  laminated  2-2 
piezocomposites  have  been  analyzed.  In  spite  of  the  fact  that  many  effective  parameters  of  the  material 
can  be  derived  by  this  simple  method,  one  key  issue  in  the  design  of  a  composite  transducer,  i.e.,  the 
^-influence  of  the  aspect  ratio  of  the  ceramic  plates  (or  the  unit  cell)  in  a  composite  on  the  performance  of 
the  transducer  cannot  be  addressed.  In  this  section,  we  will  treat  the  vibration  problem  of  a  finite 
thickness  composite  plate  under  an  external  driving  electric  field  and  situated  in  air.  From  the  analysis, 
one  can  obtain  detailed  information  on  how  the  surface  vibration  profile  changes  with  frequency  and  its 
dependence  on  the  aspect  ratio  of  the  ceramic  plate  (or  the  unit  cell  dimension),  the  possible  resonant 
modes  in  a  composite  transducer,  and  the  dependence  of  the  electromechanical  coupling  factor  on  the 
aspect  ratio  of  the  ceramic  plate,  etc.  We  will  also  show  that  as  long  as  the  thickness  resonance 
frequency  is  below  the  lateral  mode  frequency,  the  aspect  ratio  will  not  have  a  direct  effect  on  the 
vibration  uniformity  of  a  composite  near  the  thickness  resonance.  The  influence  is  on  the  bandwidth  in 
which  the  ceramic  and  polymer  vibrate  with  nearly  same  amplitude  and  phase. 

A  composite  plate  with  a  thickness  t  is  drawn  schematically  in  figure  1(b).  For  the  problem 
treated  here,  the  two  free  surfaces  of  the  composite  plate  are  electroded  with  conducting  material  and  an 
AC  electric  field  of  a  frequency  f  is  applied  between  the  two  electrodes. 

It  is  well  known  that  there  exist  no  simple  solutions  for  the  vibration  problem  of  a  finite 
thickness  plate  such  as  the  piezocomposite  treated  here.  One  of  the  most  frequently  used 
approximation  methods  is  to  expand  the  elastic  and  electric  fields  in  a  material  in  terms  of  the 
eigenfunctions  in  an  unbounded  one.21  Different  techniques  such  as  the  variational  technique  and  the 
method  of  least  squares  can  be  used  to  determine  the  expansion  coefficients.  For  the  composite  plate 
treated  here,  we  found  that  the  variational  technique  is  more  appropriate  in  treating  the  boundary 

problem  than  the  method  of  least  squares  because  of  the  large  difference  in  the  vibration  amplitudes 
between  the  two  phases  in  the  composite. 
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For  the  ceramic  phase  in  the  composite  plate,  the  elastic  displacements  u3  and  u,,  and  the 
electric  potential  <I>  are  expended  in  terms  of  the  eigenfunctions  derived: 

m  3 

u3  =  Y^LRnifni  C0S(Kix\ ) sin(/3n .r, )An 


«iC  =  'L'LRri8ni^(^ixl)cos(fin.xi)An  +  CT^sin^x,) 

n- 1  /=! 
m  3 

=  2  2  ^  cos(fe )  sin(^n.t3  )4  +  C,r3 


(18) 


Similarly,  for  the  polymer  phase  in  the  composite  plate: 

m  2  i 

<  =  22^/m  cosC^^  --))sin(j3nx3)4 

n=l  i=l  2 

m2  ^ 

8ni sin(^(x,  --))cos(yS„.r3)4  +  C<sin(/2opI(x1  --))  (  19) 

«=i  t=i  ^  2 

P  HI  p  rj 

^  =  SQ;  coshf/)^  (X]  —  — ))sin(/),,.r3  )A/;  +  Cx^ 

n= 1  ^ 

where  i^,2,  i^3,  /^j,  Rn2  and  Cn  are  determined  in  eq.  (16),  Pn,  and  /if  are  the  wave 

vector  components  of  the  nth  mode  in  the  x3-  and  x,-directions  where  the  superscripts  c  and  p  stand  for 
the  ceramic  and  polymer,  respectively,  and  are  related  to  the  mode  in  which  p  is  equal  to  zero: 


E  (O  and  -  -  co .  Rq  and  Rq  are  determined  by  eq.  (16)  in  which  S  is  set  to  zero.  The 

Ci  i  \ \  C.. .  1 


mode  with  P=0  is  generated  due  to  the  fact  that  the  velocity  of  the  electromagnetic  wave  is  much  faster 
than  that  of  the  elastic  waves  and  the  composite  plate  is  clamped  in  the  x3 -direction  (S3  =  0).  \  and  C 

are  the  coefficients  which  will  be  determined  by  the  boundary  conditions  which  are  traction  free  and 
0=±V/2  (here  exp(-jcot)  is  omitted)  at  x3=±  t/2. 


The  stresses  T3,  T5  and  electric  displacement  D,  are  expressed  as: 
r3C  =  I2^/)cosfe)cos^x3H+(e33  +  cos(h%xxx))C 


(20a) 


~£  —  V  X'T.Cr*,  r\*:~/uC. 


T5  =  1  1 7]  3  (n,/)sin(/;„>1)sin(/3/^3)An 


m  £ 

Di  =  22  AC  («,i)cos(/z^V|)cos(^nx3)An  +  (-£s33  +  D0C  cos  (/^x,))C 


(20b) 


(20c) 


T3P  =  22r33(«>0cos(^(x, -l))cos(Pnx3)An  +  T0pcos(h0pl(xl -~))C 

n=\  /=1  L  2 


( 20d ) 
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( 20e ) 


p  m  2  p  -y 

r5  =  S  I7i3(«,Osin(^.(x1--))sin(4.vOAI 

n=li=l  2 

m  2  i 

Di  =  XI^(»,Ocos(fe  --))cos(^.v,)4  -ef,C 

/-l  ^ 


(20f) 


are 


Where  ?£  («,/),  7£(n,z),  AC(».0,  Z£(/i,i),  7J3  (/z,z)  and  D3VO 

^33  (0*  2^(0,  Z)3C(0,  2j3(z),  7jf(z),  and  Z)3P(z)  in  equations  (11)  to  (14)  if  the 
replaced  by  n  and  i  for  i-th  partial  wave  of  n-th  mode  here.  And 

=  cfo'R oC>  T0p  =  cPXX,  =  e3XX- 


the  same  as 
subscript  i  there  is 


The  number  of  the  eigenfunctions,  m,  required  in  the  expansion  is  determined  by  the  accuracy 
needed  for  the  solution.  For  the  problem  treated  here,  we  found  that  it  is  adequate  to  use  eight  eigen- 
modes  in  the  expansion.  In  the  frequency  range  studied  (  fd<2  MHz*mm),  there  are  two  branches 
having  real  (3  and  the  other  branches  having  either  imaginary  or  complex  (3,  which  corresponds  to  the 
modes  confined  at  the  surface  of  x3=±t/2. 

For  the  problem  treated  here,  all  the  stress  components  in  air  are  zero  and  for  the  sake  of 
simplicity,  D  in  air  is  also  assumed  to  be  zero  since  the  dielectric  permittivity  of  the  composite  is  much 
higher  than  air.  Under  these  conditions,  the  variational  formula  takes  the  following  form: 

jv[(Tk,.k  -pii,)Su‘  +Dk'k8$*yv  +  JJ(-T,)(5»;  +(<D-  ^>)8D[]lS  =  0  (  21 ) 

where  the  integration  over  time  has  been  performed  to  take  into  account  of  the  complex  notations  for 
the  quantities  in  the  integrands.19'1  The  *  represents  the  complex  conjugate  of  the  corresponding 
quantity.  The  first  integral  is  over  the  volume  of  the  2-2  composite  plate  and  it  can  be  shown  that  it  is 
equal  to  zero  since  all  the  quantities  in  the  integrand  satisfy  eq.  (1).  The  second  integral  is  over  the 
surfaces  of  the  composite  plate  at  x3=±t  /2  where  <E>=±V/2  at  x3=±t/2  where  V  is  the  applied  voltage. 
The  periodic  condition  of  the  composite  in  the  xx  -direction  and  the  symmetric  condition  of  the  solutions 
and  the  boundary  conditions  about  the  plane  of  x3=  0  allow  the  second  integration  to  be  performed 
over  one  unit  cell  at  x3  =  t  /2. 

Substituting  T3,  T5,  upu3,  O  and  D3  into  equation  (21)  yields  the  following  linear  algebraic 
equation: 

{Mij){Aj)  =  (Vi)  (22) 

where  (My)  is  a  9x9  matrix,  (AJ)=(A1,A2,A3,A4,AJ,A6.A7,A8,C)T,  (VV)  is  a  9x1  matrix  whose  elements 
depend  only  on  V,  the  applied  voltage.  For  a  given  frequency  f,  one  can  solve  eq.  (22)  to  obtain  Aj 
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and  C.  From  A,  and  C,  eqs.  (18),  (19),  and  (20)  yield  all  the  characteristic  properties  related  to  the 
vibration  of  a  2-2  composite  plate  with  different  thickness  and  different  ceramic  volume  fraction  such 
as  the  electrical  impedance,  surface  displacement  distribution,  resonant  modes,  electromechanical 
coupling  coefficient,  etc. 


For  the  composite  plate  treated  here,  the  electrical  impedance  for  a  single  repeating  unit  can  be 
found  from: 


(23) 


where  I  is  the  current  which  is  equal  to  I  =  ^  =  -jcoQ  =  -jlcob  J  D,dxv  where  Q  and  b  are  the 


dA 


electric  charge  and  the  length  in  the  ^-direction  of  the  plate  respectively.  To  compare  with  the 

Experimental  result,  the  current  I  should  be  multiplied  by  N,  the  number  of  repeating  unit  in  a 
composite  plate. 


Shown  in  figure  6(a)  is  the  electric  impedance  curve  calculated  from  eqs.  (22)  and  (23)  for  a 
composite  plate  with  t/d=4.5  and  the  ceramic  volume  fraction  v=44%.  The  electric  impedance 
measured  experimentally  from  the  same  composite  plate  is  shown  in  figure  6(b)  and  clearly  the 
theoretical  impedance  curve  reproduces  the  experimental  data  quite  well.  The  difference  in  the 
sharpness  of  the  resonant  peaks  between  the  experimental  data  and  the  theoretical  curve  is  due  to  the 
fact  that  in  the  theoretical  analysis,  the  electrical  and  mechanical  losses  of  the  ceramic  and  polymer 
phases  were  not  included. 

In  figure  6(c),  the  peak  positions  from  the  experimental  data  are  compared  with  the  dispersion 
curves  for  this  composite  which  shows  excellent  agreement  between  the  two.  In  order  to  elucidate  the 
nature  of  these  modes,  the  spatial  distribution  of  u3  at  each  mode  at  the  composite  surface  is  presented 
in  figure  7  which  are  evaluated  based  on  eqs.  (18)  and  (19).  Apparently,  fLI  is  the  fundamental 
thickness  resonance  and  ftI  is  the  first  lateral  mode  as  revealed  by  the  fact  that  the  ceramic  and  polymer 
vibrate  180°  out  of  phase  at  this  mode,  as  predicted  in  the  earlier  theoretical  work.7’8  The  frequency 
position  and  the  distribution  of  u3  along  the  x3  axis  indicate  that  fL3  is  the  third  harmonic  of  the 
thickness  mode.  However,  the  appearance  of  ft2  is  not  expected  from  the  earlier  theoretical  works  in 
which  the  ceramic  and  polymer  vibrate  in  phase.  By  examining  the  equations  of  the  boundary 
conditions  at  x3=±t/2,  it  can  be  deduced  that  a  resonance  will  occur  whenever  3  =  (l+2n)7t/t,  i.e.,  cos 
(Pt/2)  =  0.  Hence,  the  dispersion  curves  of  real  (3.  as  shown  in  figure  6(c),  reveal  that  the  fundamental 
thickness  resonance  and  the  first  lateral  resonance  occur  at  (3  =  n/t  (fL1  and  ftl).  Similarly,  when  p  =  3 
7t/t,  the  third  harmonic  of  the  thickness  mode  will  occur  at  fL3.  In  addition,  a  mode  ft2  will  also  show 
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Figure  6.  The  electric  impedance  curve  for  a  2-2  composite  plate  with  PZT-5H  and  Spurr  epoxy 
(44%  ceramic)  measured  in  air:  (a)  the  theoretical  curve  and  (b)  the  experimental  curve.  The 
thickness  t  of  the  composite  is  t/d  =  4.5  and  d=0.635  mm.  fL1  and  fL3  are  the  fundamental  and  third 
harmonic  of  the  thickness  mode  and  tl  and  t2  are  the  modes  arose  from  the  periodicity  of  the 
composite,  (c)  The  dispersion  curves  elucidating  the  origin  of  the  resonant  modes  f„  and  ft2.  In 
general,  any  modes  on  the  dispersion  curves  will  show  up  whenever  Bd  =  (2n+l)/2  n=0  12'  is 
satisfied,  where  (3=2ji/X  and  X  =2t. 
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Figure  7.  The  distribution  of  the  elastic  displacement  u3  for  fL1,  f  L3,  f  tl,  and  f  t2  of  a  2-2  composite 
with  t/d=4,  where  (a)  and  (b)  are  the  distributions  along  the  x  3-direction  when  x  ,=0  (at  the  center 
line  of  the  ceramic  plate,  figure  1(b)),  x3=0  is  at  the  center  and  x3=2  is  at  the  ceramic  surface;  (c)  and 
(d)  are  surface  vibration  profiles  at  x3=2,  x,=0  corresponds  to  the  center  of  the  ceramic  plate  and  x 
,=0.5  is  at  the  center  of  the  polymer  plate. 
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up  at  the  branch  1  which  is  at  a  frequency  near  and  above  ftl.  By  the  same  argument,  it  would  be 
expected  that  f^,  ft3,  etc.  may  also  be  observed,  depending  on  the  electromechanical  coupling  factors 
of  these  modes.  As  shown  later,  the  effective  coupling  factor  for  the  modes  in  the  first  branch  will 
decrease  with  increasing  pd,  i.e.,  reducing  thickness,  and  on  the  other  hand,  the  effective  coupling 
factor  for  the  modes  in  the  second  branch  will  increase  with  (3.  As  the  ratio  of  t/d  decreases,  the 
frequency  of  the  thickness  mode  will  gradually  move  towards  the  first  lateral  mode  which  will  become 
stronger  (coupling  factor  increases),  and  the  second  lateral  mode  (ft2)  becomes  weaker  and  finally  it 
will  disappear.  These  features  are  summarized  in  figure  8  which  provide  understanding  on  the  earlier 
experimental  observations  on  how  various  modes  change  with  temperature  (which  causes  reduction  of 
the  shear  velocity  of  the  polymer  phase)  and  the  composite  thickness.5 


Both  the  experimental  results  and  the  theoretical  data  indicate  that  the  ceramic  and  polymer 
vibrate  in  phase  for  all  the  modes  on  the  first  branch  and  out  of  phase  for  modes  on  the  second  branch. 


The  electromechanical  coupling  factor  for  the  thickness  resonance  can  be  evaluated  based  on 
the  definition  of  IEEE:20 


k; 2  =  -£-tm (JLL-Jl) 

-  fP  2  / 


(24) 


where  k,  is  the  thickness  mode  coupling  factor,  fs  and  fp  are  the  series  and  parallel  resonance 
frequencies,  respectively.  Eq.  (24)  is  used  here  to  calculate  the  coupling  factor  for  the  modes  in  both 
the  first  and  second  branches.  Shown  in  figure  9  is  the  results  for  a  2-2  composite  plate  with  44% 
ceramic  volume  content  for  different  d/t,  where  both  theoretical  and  experimental  results  are  presented. 
As  the  ratio  of  d/t  increases,  the  coupling  factor  of  the  mode  (thickness  mode)  in  the  first  branch 
gradually  decreases,  while  the  coupling  factor  of  the  mode  in  the  second  branch  gradually  increases 
due  to  the  modes  coupling.  As  d/t  increases  further,  the  thickness  mode  will  jump  to  the  second  branch 
when  the  coupling  factor  in  the  second  branch  surpasses  that  of  the  first  branch.  Although  in  this 
region  the  coupling  factor  for  the  thickness  mode  can  still  be  quite  high,  the  distribution  of  u3  is  not 
uniform  on  the  composite  surface  and  the  ceramic  and  polymer  vibrate  180°  out  of  phase  which  is  not 
desirable  since  the  polymer  phase  will  not  be  able  to  perform  properly  the  function  of  transferring  the 
acoustic  energy  between  the  ceramic  plates  and  the  external  medium. 


For  a  composite  plate  to  work  effectively  as  an  electromechanical  transduction  material,  it  is 
required  that  the  ceramic  and  the  polymer  plates  in  the  composite  vibrate  in  phase  with  nearly  the  same 
amplitude  in  the  x3-direction.  The  evolution  of  the  vibration  pattern  in  the  two  phases  with  frequency 
and  the  effect  of  the  aspect  ratio  t/d  of  a  composite  plate  on  this  distribution  are  investigated.  Shown  in 
figure  10(a)  is  the  variation  of  the  ratio  wf  /nf  at  the  surface  of  the  composite  plate,  where  and  wf 
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Figure  8.  The  evolution  of  various  modes  in  a  2-2  composite  made  of  PZT-5H  (44%  ceramic)  and 
Spurr  epoxy  with  the  thickness  t  of  composite  plate:  (a)  t/d  =  4,  (b)  t/d  =  2,  and  (c)  t/d  =  1.  There  is 
already  substantial  coupling  between  the  thickness  and  lateral  modes  at  t/d  =  2  for  this  composite.  At 
low  t/d  values  (thin  samples),  the  lateral  mode  will  disappear. 
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Figure  9.  The  evolution  of  the  coupling  factor  for  the  modes  on  the  first  and  second  branches  of  the 
dispersion  curves  with  the  composite  thickness  d/t  for  a  composite  made  of  PZT-5H  and  Spurr  epoxy 
with  44%  ceramic  content. 
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Figure  10.  (a)  The  ratio  of  u3p/u3c  vs.  frequency  for  a  composite  made  of  PZT-5H  and  Spurr  epoxy 
(44%  ceramic  content)  at  a  thickness  of  t/d  =  4  (d=0.635  mm)  measured  in  air.  u3p  and  u3c  are  the 
surface  displacements  at  the  centers  of  the  polymer  plate  and  ceramic  plate,  respectively.  The  black 

dots  are  the  experimental  data  measured  using  a  laser  dilatometer  and  the  solid  line  is  derived  from 
the  model. 

(b)  The  ratio  f,/fs  as  a  function  of  the  ratio  d/t  of  the  composite  plate,  where  f,  is  the  frequency  at 
which  u3p/u3c=l  and  fs  is  the  series  resonant  frequency  for  the  thickness  mode.  The  black  dots  are 
experimental  data  and  the  solid  lines  are  derived  from  the  model  calculation. 

(c)  The  bandwidth  as  a  function  of  the  ratio  d/t  of  the  composite  plate.  The  black  dots  are 
experimental  data  and  the  solid  tines  are  derived  from  the  model  calculation. 
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are  u3  at  the  centers  of  the  polymer  (x,=d/2)  and  the  ceramic  plates  (x,=0),  respectively,  with  frequency 
for  the  composite  plate  of  t/d=4  where  both  the  experimental  results  and  theoretically  calculated  curve 
are  presented.  The  experimental  data  were  acquired  using  a  laser  dilatometer.24  Hence,  at  frequencies 
far  below  any  resonant  mode,  u\  /  w3c  is  always  less  than  one,  which  is  true  as  long  as  the  composite  is 
driven  electrically.  For  a  given  frequency,  as  t/d  increases,  this  ratio  increases  and  approaches  one. 
These  are  consistent  with  the  results  of  the  earlier  theoretical  work  on  the  static  properties  of 
composites.  As  frequency  increases,  the  ratio  / zif  increases  towards  one.  At  a  frequency  f, 
which  is  near  fs  of  the  thickness  mode,  m3p/u3c=  1.  This  ratio  will  surpass  one  as  the  frequency  is 
further  increased.  This  is  true  as  long  as  fL1<ftI.  In  figure  10(b),  the  change  of  f/fs  vs.  the  ratio  of  d/t 
is  presented.  Clearly,  fj/fs  is  close  to  but  larger  than  one  except  for  composite  plates  with  a  small 
aspect  ratio.  Hence,  the  aspect  ratio  t/d  does  not  have  a  significant  effect  on  the  ratio  of  I  iu  at 
♦frequencies  very  near  fs  of  the  thickness  mode,  where  iu  /  «3C  is  always  near  one  if  the  thickness  mode 
is  below  the  first  lateral  mode  frequency.  However,  it  will  affect  the  bandwidth  in  which  it?  /u%  is  near 
one.  For  example,  the  bandwidth  Af/fp  where  Af  is  defined  here  as  the  frequency  width  in  which  0.9 
<lh  /m3  <1.1,  increases  as  the  ratio  of  d/t  decreases,  which  is  shown  in  figure  10(b).  The  experimental 
data  points  are  also  presented  in  figure  10(b)  and  the  agreement  between  the  two  is  quite  good. 

2 A-  Forced  Vibration  of  a  Piezocomposite  Plate  in  a  Fluid  Medium: 

Following  a  similar  procedure  as  outlined  in  the  preceding  section,  the  vibration  problem  of  a 
composite  plate  in  a  fluid  medium  can  also  be  analyzed  which  are  more  relevant  to  the  practical  design 
and  application  of  a  composite  transducer.  The  fluid  medium  chosen  for  the  study  is  water  and  two 
situations  will  be  investigated:  the  composite  plate  as  a  transmitter,  i.e.,  under  a  harmonic  electric  field, 
and  the  composite  as  a  receiver,  i.e.,  under  a  harmonic  acoustic  pressure  in  water. 

(2.4.1)  Forced  Vibration  of  a  Composite  Plate  in  Water  under  an  AC  Electric  Field 

For  the  sake  of  convenience,  we  will  treat  the  system  as  a  symmetric  one  in  which  the 

composite  plate  is  loaded  by  water  on  both  surfaces,  hence,  the  boundary  conditions  at  the  composite- 
water  interface  are: 

TCm  _  TW  rpCm  _  n  Cm  W  ,  _  ,  V  t 

*3  ~  *3  »  -*5  ~  0,  M3  =  and  <I>  =  ±—  at  X3  =  ±—  ( 25 ) 

where  superscripts  Cm  and  W  denote  the  quantities  in  the  composite  and  water,  and  the  factor  e~ja*  is 
omitted  in  electric  potential  O.  The  variational  formula  for  this  problem  can  be  derived  as 
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(26) 


Jj-  ~[(Tf -  T3CmX6(uf)'  +  5(«fV)  +  (uf"'  -  u^xsa^f  +  d(T3Cmf)j/S 
+t[(4,c“  -W(D$m)‘  -  rf"‘SO,f”‘)'  ]ds  =  0 

where  S  indicates  that  the  integration  is  over  the  two  composite-water  interfaces  (at  x3=  ±  t/2),  and 
the  integration  over  the  volume  has  been  omitted  since  it  is  equal  to  zero. 


Here,  the  expressions  for  the  elastic  displacement  vectors  and  the  electrical  potential  O  in  the 
composite  are  those  in  eq.  (18)  and  (19).  Because  of  the  periodic  nature  of  the  composite  plate  in  the 
Xj-direction,  the  solutions  in  water  have  the  form: 

UY  =J  iC sin(hn%)exp(±y^ x3)R^ 

n= 0 

uf  =  I  cos(h^,  x{ )  exp(±jj3^x3)R^V  (27) 

n~  0 

*3  -itf  Xcos<il)exp(±j^3)i?„'v 

n= 0 


where 


h™  =  Bw  = 

d  Pn 
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IV \2 


\(vw)2 


-(K) 
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=  c’W)2 


Bw  =  — 

HQ  y\V  5 


V 


w 


is  the 


longitudinal  wave  velocity  of  water,  cw  is  the  bulk  modulus  of  water.  ±  correspond  to  the  solutions  in 
x3  >0  region  and  x3  <  0  region,  respectively,  h  and  (3  are  the  wave  vector  components  in  the  xr  and 
x3-directions  in  water.  R™  is  determined  by  the  boundary  conditions  through  the  variational  principle 
formula  eq.  (26). 


Substituting  equations  (18),  (19)  and  (27)  into  equation  (26)  yields  a  set  of  linear  algebraic 
equations: 

(Mjj)(Aj)=(Vi)  (28) 

where  (My)  is  a  matrix  in  which  its  elements  are  related  to  the  parameters  of  water,  ceramic  and 
potymer  and  the  geometrical  parameters  v  and  d,  as  well  as  co,  (3  and  h.  (Aj)=(AI,  ...,  A  ,  C, 
^  w^ere  m  and  J  are  the  numbers  of  the  eigenfunctions  used  in  the  expansions  for  the 
quantities  in  the  composite  and  in  water,  respectively,  and  how  many  eigenfunctions  should  be  used  in 
the  expansion  depends  on  the  accuracy  desired.  In  this  calculation,  m=8  and  J=6.  Hence,  (V,)  is  a 
16x1  matrix  which  elements  depend  on  the  applied  voltage  V.  Equation  (28)  is  solved  numerically. 


Shown  in  figure  11(a)  is  a  theoretical  electrical  impedance  curve  for  a  2-2  composite  plate  of 
ceramic  volume  content  and  the  aspect  ratio  t/d  =  3.8  loaded  with  water.  For  the  comparison,  the 
xperimental  curve  for  the  same  composite  plate  is  shown  in  figure  1 1(b).  The  agreement  between  the 
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experimental  result  and  the  theoretical  one  is  quite  good.  The  relatively  sharp  resonant  peaks  in  the 
theoretical  curve  compared  with  the  experimental  one  are  due  to  the  fact  that  in  the  theoretical 
calculation,  the  elastic  and  dielectric  losses  of  the  composite  plate  are  not  included.  Comparison 
between  figure  1 1  with  figure  6,  which  is  the  electrical  impedance  for  a  similar  composite  plate  in  air, 
reveals  that  the  resonance  is  severely  damped  in  water  as  shown  by  the  marked  broadening  of  the 
resonant  peaks  in  the  impedance  curve.  For  the  lateral  modes  (modes  1  and  2  in  figure  11)  which  is 
mainly  determined  by  the  shear  resonance  of  the  polymer  phase  in  the  2-2  composite  investigated, 
because  the  acoustic  impedance  between  the  polymer  and  water  is  very  close,  the  change  in  the 
amplitude  is  quite  significant.  One  interesting  feature  of  the  influence  of  the  water  loading  on  a 
composite  plate  is  the  appearance  of  a  mode  which  is  labeled  3  in  the  figure.  It  can  be  shown  that  the 
resonant  frequency  for  this  mode  does  not  change  very  much  as  the  thickness  of  the  composite  plate 
^changes.  From  the  stress  distribution  pattern  in  the  water,  it  is  not  difficult  to  show  that  this  mode  is 
related  to  the  coupling  between  the  composite  and  water  and  its  frequency  is  determined  by  the 
periodicity  d  and  the  acoustic  wave  velocity  of  water  (f  =  V7d).  And  it  corresponds  to  the  local 
oscillation  of  water  within  one  unit  cell. 

In  analogous  to  the  situation  in  air,  it  is  also  found  that  in  water,  at  frequencies  far  below  any 
resonant  mode,  the  ratio  of  w3P  /«3C  is  always  less  than  one  and  will  approach  one  as  the  aspect  ratio  t/d 
increases.  As  the  frequency  increases  towards  the  thickness  resonant  frequency,  w3P/n3c  increases 
towards  one.  At  a  frequency  f,  near  the  thickness  resonant  frequency  fs,  this  ratio  becomes  one  and 
above  that,  this  ratio  is  larger  than  one  (figure  12(a)).  Figure  12(b)  presents  f/fs  and  Af/f,  vs.  the  ratio 
of  d/t  for  a  composite  with  44%  ceramic  content.  Clearly,  the  effect  of  the  aspect  ratio  t/d  on  the 
surface  uniformity  of  a  composite  plate  is  to  change  the  frequency  position  f,  with  respect  to  the 
thickness  resonance  and  the  frequency  width  Af  in  which  lu  /uf  is  near  one  which  are  very  similar  to 
that  found  in  air.  However,  the  bandwidth  Af/f,  in  water  is  larger  than  that  in  air,  indicating  that  the 
surface  distribution  of  the  displacement  of  the  two  phases  is  much  flatter  in  water  than  in  air  which  is 
quite  understandable.  Because  of  the  water  loading,  the  vibration  amplitude  of  the  polymer  phase  is 
significantly  reduced.  Figure  12(a)  shows  how  the  ratio  of  u3p/u3c  at  the  composite  surface  varies  with 
frequency  (t/d=4).  One  noticeable  change  between  the  surface  profiles  in  air  and  in  water  is  that  in 
water,  even  at  the  lateral  mode  frequency,  this  ratio  does  not  become  very  large  (not  shown  in  the 
figure  since  the  lateral  mode  frequency  is  at  about  1.5  MHz). 


26 


Figure  11.  The  electrical  impedance  magnitude  for  a  2-2  composite  of  PZT-5H  and  Spurr  epoxy  with 
44%  ceramic  content  measured  in  water.  The  thickness  of  the  composite  is  t/d=3.8  and  d=0.635 
mm.  (a)  is  derived  from  the  model  and  (b)  is  the  experimental  results.  The  modes  1  and  2  are  the 
lateral  modes  (f„  and  ft2  in  figure  6(c))  and  mode  3  arises  from  the  coupling  of  water  to  the  periodic 
structure  of  the  composite  surface. 
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Figure  12.  (a)  The  ratio  of  u3p/u3c  vs.  frequency  for  a  composite  made  of  PZT-5H  and  Spurr  epoxy 
(44%  ceramic  content)  at  a  thickness  of  t/d  =  4  (d=0.635  mm)  measured  in  water.  u3p  and  u3c  are  the 
surface  displacements  at  the  centers  of  the  polymer  plate  and  ceramic  plate,  respectively.  The  results 
are  derived  from  the  model. 

(b)  The  ratio  f,/fs  as  a  function  of  d/t  of  the  composite  plate,  where  f,  is  the  frequency  at  which 
u3p/u3c=l  and  fs  is  the  series  resonant  frequency  for  the  thickness  mode.  The  results  are  derived  from 
the  model. 

(c)  The  bandwidth  as  a  function  of  d/t  of  the  composite  plate.  The  results  here  can  be  compared 
with  those  in  figure  10  and  apparently,  and  the  water  loading  improves  the  uniformity  of  the 
vibration  profile  at  the  composite  surface.  The  results  are  derived  from  the  model. 


28 


(2.4.2)  Forced  Vibration  of  a  Composite  Plate  in  Water  under  Harmonic  Acoustic  Pressure 


Now  we  turn  to  investigate  the  vibration  behavior  of  a  composite  plate  in  water  under  harmonic 
acoustic  pressure,  that  is,  to  study  issues  related  to  the  receiving  sensitivity  of  a  simple  2-2  composite 
transducer.  For  simplicity,  we  shall  restrict  the  treatment  to  the  symmetric  system.  In  this  case,  as 
sketched  in  figure  13,  the  incident  acoustic  waves  impinge  normally  on  the  two  surfaces  of  the 
composite  plate  from  opposite  directions  and  the  boundary  conditions  for  this  problem  are: 


T-fm  =  7f ,  Tfm  =  0,  u%m  =  uf  and  D%m  =  0  atx3=±|  (29) 


l3  ~  v3  -  J5 
The  appropriate  variational  formula  is 


t  |[(3f  -  Tfmmu¥?  +  S(ufmf)  +  (ufm  -  uf  )(6(T3w )*  +  5(Tfmffys 
-Jj.  [E$mS(<bCmf  +  TfmS(uf"'y  ]dS  =  0  ( 


The  expressions  of  the  displacements  and  electrical  potential  in  the  composite  plate  are  those  in 
equations  (18)  and  (19)  except  now  the  mode  of  [3=0  should  not  be  included  since  there  is  no 
externally  applied  electric  field.  To  account  for  the  incident  wave  in  water,  equation  (27)  is  modified 
and  the  solutions  in  the  x3  >  0  region  are: 

UT  =  J  lhn  sin (C *i )exp (jp™x3 )R™ 

n=  0 

uf  =  -p0  exp HP0X3 )  +  HPn  CO s(h^xl )  exp(y'/?^ x3  )R^  (31) 

n= 0 

=  JT0V [exp(— y'/?0-v3 )  +  X  cos(/z7JV.v1 )  exp( jP^x3  )R™  ] 
n= 0 

where  the  term  of  ~P0  exp(-jpQx3 )  corresponds  to  the  incident  plane  wave.  By  substituting 

equations  (18),  (19)  and  (31)  into  (30),  we  obtain  a  set  of  linear  algebraic  equations* 

(Mij)(Aj)=(Vi)  ( 32 ) 


where  (My)  is  a  15x15  matrix,  (Aj)=(A1,...,  A^,  RqV ,  ...,  Rw )  with  m=8  and  J=6,  and  (Vt)  is  a  15x1 
matrix  related  to  the  incident  wave.  By  solving  equation  (32)  one  can  obtain  all  of  the  properties  related 
to  the  behaviors  of  a  composite  plate  under  a  harmonic  acoustic  pressure. 


Shown  in  figure  14  is  the  open  circuit  voltage  receiving  sensitivity  V/p,  where  V  is  the  voltage 
output  and  p  is  the  pressure  of  the  incident  wave,  for  a  44%  2-2  composite  plate  at  the  thickness 
resonance  vs.  d/t  for  the  thickness  resonance  mode  where  the  cross  over  region  (point  B  in  figure  2)  is 
at  d/t  —  0.65.  That  is,  at  d/t  <  0.65,  the  thickness  mode  is  in  the  first  branch  of  the  dispersion  curves 
and  above  that,  the  thickness  resonance  is  at  the  second  branch  of  the  dispersion  curves  where  the 
polymer  and  ceramic  vibrate  out  of  phase.  Since  the  problems  treated  in  this  paper  are  related  only  to 
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the  piezo-materials,  the  open-circuit  receiving  sensitivity  is  used  here  which,  as  pointed  out  by  Kojima, 
is  a  system  independent  parameter.26  The  bandwidth  is  defined  as  the  3dB  width  about  the  peak 
frequency  (fp  for  the  receiving  mode).  For  a  single  phase  material,  it  can  be  derived  from  the  KLM 
model  that  the  sensitivity  here  should  be  proportional  to  t,  the  thickness  of  the  transducer.  Hence,  in 
figure  14,  the  V/p  vs.  d/t  curve  for  a  single  phase  material  should  fall  off  as  t  (solid  line  in  the  figure) 
while  the  bandwidth  should  stay  constant.  The  results  show  that  the  sensitivity  of  the  thickness  mode 
for  a  2-2  composite  decreases  slowly  as  d/t  becomes  larger  than  0.4,  but  the  bandwidth  increases 
gradually.  The  increase  in  the  bandwidth  is  due  to  the  merger  of  the  two  resonant  modes  in  water. 
After  that,  there  are  anomalous  changes  in  both  the  bandwidth  and  the  sensitivity  in  the  cross-over 
region.  At  higher  values  of  d/t  (thin  composite  plates),  both  the  sensitivity  and  the  bandwidth  fall 
much  below  the  values  of  single  phase  material  (solid  line).  If  we  define  the  figure  of  merit  (FOM) 
r^ere  as  the  product  of  sensitivity  and  bandwidth,  as  shown  in  figure  14(b),  at  d/t  <  0.5,  the  FOM  falls 
off  with  t.  At  d/t  >  0.5,  the  FOM  drops  to  much  smaller  value.  Therefore,  in  order  to  gain  a  high 
receiving  sensitivity  and  a  broad  bandwidth  of  a  2-2  composite  transducer,  it  is  desirable  to  have  d/t 
less  than  0.5.  The  results  here  can  be  compared  with  what  shown  in  figures  2,  8,  9,  and  12.  At  d/t 
above  0.5,  the  coupling  factor  shows  a  precipitous  drop  (figure  9)  and  the  thickness  mode  frequency 
also  shows  an  apparent  deviation  from  the  extrapolated  value  (the  dashed  line  in  figure  2).  Figure  8 
also  shows  that  there  is  a  significant  coupling  between  the  two  modes.  Although  the  results  presented 

are  for  a  composite  with  44%  ceramic  content,  it  is  approximately  true  for  composites  with  other 
volume  fractions. 

In  general,  the  maximum  transmitting  voltage  sensitivity  of  a  transducer  is  at  a  frequency  near 
fs  and  the  maximum  open  circuit  receiving  sensitivity  is  near  fp  of  the  thickness  mode  (fp  >  fs),  and 
hence  it  is  expected  that  the  influence  of  the  lateral  mode  (hence,  the  aspect  ratio  t/d)  will  be  more 
severe  on  the  receiving  sensitivity.  From  the  data  analysis,  it  is  shown  that  when  an  FOM  which  is  the 
product  of  the  sensitivity  and  the  bandwidth  is  introduced  as  the  criterion  of  the  performance  of  a 
composite  transducer  as  a  receiver,  the  performance  deteriorates  when  d/t  >  0.5,  for  the  composite 
discussed  here.  Experimental  results  confirm  this  finding  where  as  shown  in  figure  15,  the  insertion 
loss  of  a  2-2  piezocomposite  shows  large  increases  at  d/t  near  and  above  0.5.27 
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Figure  13.  Schematic  drawing  of  a  2-2  composite  plate  under  normal  incident  of  a  pressure  wave 
from  the  water. 


Figure  14.  (a)  The  influence  of  the  unit  cell  dimension  ratio  of  d/t  on  the  open  circuit  voltage 
sensitivity  of  a  2-2  composite  with  44%  PZT-5H  and  Spurr  epoxy  matrix.  Solid  dots  are  V/p  at  the 
peak  frequency  (fp)  and  the  solid  line  is  the  sensitivity  of  a  single  phase  material  with  the  effective 
properties  of  the  composite.  The  results  are  derived  from  the  model. 


(b)  The  3  dB  bandwidth  as  a  function  of  d/t  for  a  composite  transducer  derived  from  the  model. 

(c)  FOM  (solid  dots,  the  product  of  the  sensitivity  and  the  bandwidth)  for  the  2-2  composite  in  (a) 
as  a  function  of  d/t.  The  solid  line  is  FOM  for  a  single  phase  material  with  the  effective  properties  of 

the  composite.  At  d/t  >  0.5,  FOM  drops  much  below  the  single  phase  material  value.  The  results 
are  derived  from  the  model. 


31 


*  Jhe  exPenmental  data  on  insertion  loss  (open  circles)  and  6  dB  bandwidth  (solid 

fA2  Cu°mp?lteS  Wlth  44%  volume  content  of  PZT-5H  and  Spurr  epoxy  polymer  matrix 

d^wn?d  by^the  PUl$e  £Ch0  meth°d  f°r  different  ratio  of  W-  (d=0.635  mm).  The  solid  lines  are 
drawn  to  guide  eyes. 
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HI  Plane  Wave  Propagation  at  an  Interface  of  a  Medium-2-2  Piezocomnosife  and  the  FfferrivP 
Impedance  of  the  Composite28,29 


In  this  section,  the  piezoceramic  polymer  composites  used  in  the  model  calculations  and  in 
experiments  are  those  made  of  PZT-5H  piezoceramic  and  Spurrs  epoxy  polymer  matrix  (see  Table  I). 


3-1-  Pl<ms  acoustic  wave  propagation  at  a  fluid-2-2  composite  interface: 


(3.1.1)  Eigen-modes  and  eigen-functions  in  the  2-2  composite  and  fluid  medium 


When  a  plane  acoustic  wave  encounters  a  boundary  between  two  media,  energy  is  reflected 
and  transmitted  from  and  across  the  boundary.  In  this  section,  we  shall  derive  the  former  formula 
used  to  treat  the  plane  wave  transmission  and  reflection  at  the  interface  between  a  fluid  and  a 
piezoceramic  polymer  composite  of  2-2  connectivity.  The  schematic  of  the  interface  system  is 
shown  m  figure  16  where  a  plane  acoustic  wave  is  normally  incident  upon  the  boundary  between  a 
fluid  and  a  2-2  piezocomposite.  The  coordinate  system  is  chosen  such  that  the  x3-axis  is 
perpendicular  to  the  interface,  the  Xj-axis  is  parallel  to  the  interface  as  shown  in  the  figure,  and  the 
x2-axis  is  perpendicular  to  both  the  x,  and  x3  axes  (pointing  out  of  the  paper).  The  fluid  occupies 
the  space  of  x3  <  0  and  the  2-2  composite  is  in  x3  >  0.  To  the  problems  treated,  apparently,  the 
periodic  lamina  structure  of  the  composite  in  the  xrdirection  causes  variations  of  the  acoustic  and 


incident  plane  wave 


laminated 

composite 

ceramic 

polymer 


Figure  16.  Schematic  of  an  interface  between  a  fluid  and  a  2-2  piezocomposite  where  a  plane 
acoustic  wave  normally  incidents  to  the  interface  from  the  fluid  medium.  The  composite 
occupies  the  semi-infinite  space  x3  >  0. 
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electric  fields  in  that  direction  and  the  interface  between  the  composite  and  fluid  introduces 
variations  of  these  fields  in  the  x3-direction. 


The  eigenmodes  and  eigen-functions  for  a  2-2  piezocomposite  have  been  derived  in  the 
section  2.1  and  section  2.2. 


The  boundary  conditions  at  the  fluid-2-2  composite  interface  for  the  problem  are 
T3C3m  =  Tg,Tlc3"'=0,u$m=u? 

®Cm  =  EfCm  =  jjW  at  x3-°  (33) 


The  superscripts  Cm  and  W  represent  the  2-2  composite  and  the  fluid.  And  <t>  is  the  electric  potential. 


In  the  fluid  medium,  the  constitutive  relation  has  the  form 

r1,  =  r33=A|!.+|i) 

oX\  ox 3 


(34) 


where  A  is  the  bulk  modulus  of  the  fluid.  The  periodic  condition  in  the  x, -direction  due  to  the 
composite  structure  also  imposes  the  constraints  on  the  solutions  in  the  fluid  region- 
u1w=jhnwsin(hnwxI)exp(jP^x3)Rnw 

<  =  pnw  cos(hnwx,)exp<jP»xi)R»  (35) 


2xiTC  (^2 

where  hn  ~  anc^Pn  -  J /yWx2  -  (^n  )”  (n=0, 1,2,3...)  are  the  wave  vector  components  in  the 


d  "*rn  l|(Vw)2 


iW 


xr  and  x3 -directions  m  the  fluid  region,  Vw  =  y—  is  the  longitudinal  wave  velocity  of  fluid,  and 

the  superscript  W  denotes  the  quantities  in  the  fluid  medium.  In  eq.  (35),  Rnw  is  an  undetermined 
constant. 


(3.1.2)  Variational  formula  on  the  boundary  problem  at  the  composite-fluid  interface 


In  order  to  satisfy  the  boundary  conditions  at  the  composite-fluid  interface,  the  elastic  and 
electric  fields  in  the  two  regions  are  expanded  in  terms  of  the  eigenfunctions  obtained.  The  coefficient 

of  each  mode  is  determined  by  the  variational  method  as  briefly  outlined  in  this  section.  The  variational 
formula  for  the  problem  here  is: 

+  8u^m )*  -  2 

+  (upn  - uf  xmff  +  8rg?  +{D^  -D%mX8®w  +  8®Cmf  (36) 

+  (®W  -  0Cm)(SD f  +  SD%m)*]dS  =  0 
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d  d 

where  the  surface  integral  is  performed  at  x3=0  over  one  unit  cell  (-- -  <  x1  <  — ),  W  and  Cm  refer  to 

2  2 

the  fluid  and  composite,  respectively,  and  *  denotes  the  complex  conjugate. 

For  the  situation  considered  here,  in  the  fluid  medium,  there  are  incident  wave  and  reflected 
wave,  therefore,  in  terms  of  the  eigenfunctions  in  the  fluid  medium,  the  elastic  displacement  and 
electric  potential  in  the  fluid  can  be  expressed  as 

**l  =jl  hn  sin (h*. ^)exp (-jff. x3 )R™ 
n= 0 

q  (37) 

u3  =  Po  zxv(jPo'x3)-  X  pY  cos(h?Xl ) exp(-y'/^Vx3 )rW 

n=0 

t  =  1  cos(/z,f  Xj )  expf/z^  x3)R™ 

n= 0 

where  and  0q  exp(jpQX3)  corresponds  to  the  incident  wave  normal  to  the  interface. 

are  the  coefficients  determined  by  equation  (36). 


In  the  composite  region,  only  the  transmitted  wave  exists.  However,  the  expressions  for  the 

elastic  and  electric  variables  in  the  ceramic  phase  and  in  the  polymer  phase  are  different.  In  the  ceramic 

phase,  the  elastic  displacements  and  the  electric  potential  are 
(*  Ttl  3  ^  ~ 

“3=1  I  Rnifni  cosihpiXy )  exp(jpnx3  )An 
n=li'=l 

c  m  3  r  r  r 

“1  =j  I  1  RniSni  sin(/4'.ti )  exp{jpnx3  )An 

r  m  3  r  r 

^  ~~  X  X  RniSni  cos(fyzr^l )  CXp (7/^X3  )An 

n=li=l 

and  in  the  polymer  phase,  they  are 

p  m  2  p  p  p  si 

“3  =  I  IKifni  cos(hni(xl--))exp(jPnX3)An 


n=lz=l 

U1  =jl  I  Rnigfu  sin(^  (*1  -  77))  exp(jpnx3  )An 

n=li=l 


m 


°  “  I  cn  cosh (pn  (xx  -  -))  exp UPnx3  )An 

n-\  ± 


(38) 


(39) 


where  A„  are  coefficients  to  be  determined  from  eq.  (36).  In  the  expansions  (37),  (38)  and  (39),  the 
number  of  modes  used  are  determined  by  the  accuracy  needed.  Making  use  of  the  constitutive 
equations  in  the  ceramic  and  polymer  phases,  the  stress  and  electric  displacement  in  the  composite 
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region  can  be  obtained.  Substituting  these  quantities  into  eq.  (36)  yields  the  following  algebraic 
equation: 

(My.)Rj)  =  (Vi)  (40) 

where  (Rj)-(R0,  R,,...RG,  A,,...Am)T,  (My)  is  a  (G+m+l)x(G+m+l)  matrix,  and  (Vp  is  a  (G+m+l)xl 

matnx.  In  the  problem  treated,  G=7  and  m=9  are  used.  From  eq.  (40),  the  coefficients  R.  and  A,  are 
determined.  J  0 


(3.1.3)  Reflection  coefficient  and  input  acoustic  impedance 

For  an  acoustic  transducer  to  function  effectively,  it  is  highly  desirable  that  nearly  all  of  the 
^acoustic  energy  aniving  at  the  interface  be  transmitted  across  the  boundary,  which  implies  a  minimum 
m  the  reflected  wave.  In  this  section,  we  will  first  derive  the  reflection  coefficient  for  a  composite  in 
contact  with  water  based  on  the  results  derived  in  the  preceding  sections.  As  we  have  pointed  out,  in 
the  expression  (37),  pff  expCy/Jo x3)is  the  incident  wave  and  other  terms  are  from  the  reflected  wave. 
The  zeroth  order  reflected  wave,  i.e.,  n=0  and  h0w=0,  represents  the  normal  reflection.  For  n*0,  in 
the  frequency  of  mterest,  i.e.,  below  the  first  lateral  mode  frequency,  it  can  be  shown  that  pnw  is  nearly 
equal  to  j  (2mt/d)  which  implies  that  the  amplitude  of  these  waves  will  decay  as  exp(-2mtx3/d)  as  the 
waves  depart  from  the  interface  and  hence  they  are  evanescent  waves  which  are  important  only  in  the 

boundary  region.  In  this  sense,  the  reflection  coefficient  which  can  be  measured  experimentally  is 
equal  to  -Rq. 

Figure  17(a)  shows  the  reflection  coefficient  derived  for  a  plane  acoustic  wave  on  the  boundary 
between  water  and  a  2-2  composite  in  which  the  volume  fraction  of  ceramic  phase  is  44%.  At  lower 
frequencies  where  the  periodicity  d  of  the  2-2  composite  in  the  x,-direction  is  much  smaller  than  the 
wavelength  of  the  incident  wave,  the  2-2  composite  behaves  like  an  effective  homogeneous  medium 
and  the  reflection  coefficient  is  real  and  can  be  fully  characterized  by  the  acoustic  impedance  of  the  two 
materials,  i.e.,  the  effective  acoustic  impedance  of  the  composite  derived  based  on  effective  medium 
theory  and  the  characteristic  impedance  of  water.”0  As  the  frequency  increases,  however  the 
vibrational  displacement  at  the  surface  of  the  2-2  composite  becomes  nonunifoim,  which  implies  that  it 
could  no  longer  be  regarded  as  a  homogeneous  medium.  Hence,  the  reflection  coefficient  at  the  water 
and  composite  interface  will  be  dispersive  and  become  a  complex  number  as  shown  in  figure  17(a) 
w  ere  the  magnitude  of  the  reflection  coefficient  decreases  gradually  and  meanwhile,  the  phase  of  the 
reflection  coefficient  becomes  larger  than  180°  and  increases  with  frequency.  To  verify  the  result,  an 
experiment  was  conducted  to  measure  the  reflection  coefficient  from  the  interface  between  water  and  a 
2-2  composite  of  44%  ceramic.  The  schematic  of  the  experimental  set-up  is  shown  in  figure  18,  where 
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Figure  17.  The  reflection  coefficient  (both  the  amplitude  and  phase)  as  a  function  of  frequency  for  2- 
2  composites  with  (a)  44%  ceramic  and  (b)  15%  ceramic  content.  The  anomalous  changes  at  fd  =  1 
in  (a)  and  fd  =  0.8  in  (b)  are  related  to  the  lateral  resonant  mode  in  these  composites,  respectively. 
The  experimental  data  are  shown  in  (a)  where  the  open  circles  are  the  amplitude  and  solid  circles  are 
the  phase  of  the  measured  reflection  coefficient. 
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Figure  18.  Schematic  of  the  set-up  used  to  measure  the  reflection  coefficient  of  2-2  composites  In 
order  to  obtain  phase  information  of  the  reflection  coefficient,  the  distance  between  the  probe 

transducer  and  the  2-2  composite  is  maintained  the  same  as  that  between  the  probe  and  the  stainless 
steel  standard. 
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the  reflected  signals,  both  the  amplitude  and  time  delay,  from  a  2-2  composite  of  44%  ceramic  volume 
and  a  standard  reflector  made  of  a  stainless  steel  rod,  were  compared.  Since  the  acoustic  impedance  of 
the  stainless  steel  standard  is  know,  the  reflection  coefficient  from  the  water-composite  interface  can  be 
determined.  To  approximate  the  plane  wave  condition,  the  distance  between  the  transducer  and  the 
water-composite  interface  is  smaller  than  the  radius  of  the  transducer.  The  result  thus  obtained  is 
shown  in  figure  17(a).  The  agreement  between  the  model  result  and  experimental  result  is  quite  good. 

In  figure  17(b),  the  reflection  coefficient  for  a  2-2  composite  with  15%  ceramic  content  is 
presented.  Clearly,  the  trend  is  quite  similar  to  that  in  figure  17(a),  although  the  magnitude  and  the 
phase  of  the  reflection  coefficient  are  more  sensitive  to  the  frequency  of  incident  wave  which  is 
apparently  due  to  the  small  volume  content  of  the  ceramic  phase  in  the  composite  and  lower  lateral 
/node  frequency  (fd  ~  0.8  Mhz*mm). 

One  important  parameter  related  closely  to  the  reflection  and  transmission  of  an  acoustic  wave 
at  an  interface  is  the  acoustic  impedance  of  the  media  at  the  two  sides  of  the  boundary  which  also 
depends  on  the  nature  of  the  wave  such  as  plane  wave  or  spherical  wave.2,13  For  homogeneous 
isotropic  materials,  the  acoustic  impedance  for  a  plane  wave  is  the  characteristic  acoustic  impedance 
which  is  simply  the  product  of  the  mass  density  p  and  the  wave  (longitudinal  wave  or  shear  wave) 
velocity  V.19,31,32  For  nonuniform  materials  (such  as  piezocomposite  materials),  on  the  other  hand,  no 
characteristic  impedance  can  be  defined  in  such  a  simple  manner  due  to  the  dispersive  nature  of  the 
properties  as  has  been  shown.  However,  the  acoustic  impedance  of  a  material  can  still  be  found,  for 
example,  from  the  reflection  coefficient  or  other  methods.  It  should  be  pointed  out  that  although  there 
exist  several  definitions  for  determining  the  acoustic  impedance  at  an  interface  (input  acoustic 
impedance),  and  for  a  heterogeneous  material,  the  results  obtained  by  using  different  definitions  may 
not  be  the  same.33'34  However,  since  the  input  acoustic  impedance  is  not  a  direct  physical  quantity  but 
rather  a  parameter  introduced  for  the  convenience  of  the  analysis,  a  situation  in  analogy  to  the  complex 
notation  introduced  in  many  engineering  fields,  which  definition  is  more  appropriate  really  depends  on 
the  situation  where  it  is  used.  In  dealing  with  the  energy  transfer  across  an  interface,  one  might  have 
to  use  the  formula  derived  under  energy  flow  consideration.  On  the  other  hand,  in  designing  matching 
layers  and  dealing  with  the  reflection  and  transmission  of  acoustic  waves  in  a  multilayer  medium,  the 
amplitude  and  phase  of  the  wave  components  are  crucial  which  may  not  be  included  in  the  coefficients 
related  to  the  transmission  and  reflection  of  the  acoustic  energy.  In  the  discussion  here,  we  will  adopt 
the  definition  related  to  the  wave  reflection  and  transmission  problems.30,33 


From  the  reflection  coefficient  R,  the  input  acoustic  impedance  of  the  composite  at  the  interface 
can  be  found  as 
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(41) 


where  Zin  is  the  effective  input  impedance  of  the  2-2  composite,  Zw  is  the  characteristic  impedance  of 
water.  Using  the  results  in  figure  17,  the  effective  input  acoustic  impedance  of  2-2  composites  with 
44%  and  15%  ceramic  volume  content  is  evaluated  and  presented  in  figures  19(a)  and  19(b). 
Apparently,  there  is  a  large  change  of  the  magnitude  and  phase  of  the  input  acoustic  impedance  as  the 
frequency  increases.  For  the  comparison,  the  acoustic  impedance  is  also  evaluated  from  eq.  (41) 
where  only  the  magnitude  of  the  reflection  coefficient  is  used  and  shown  in  figure  19.  Clearly,  the 
results  show  how  erroneous  it  can  be  if  the  phase  information  in  the  reflection  coefficient  is  ignored.35 


In  analogy  to  the  electrical  impedance,  the  specific  impedance  at  the  interface  can  also  be  found 
^from  the  ratio  of  the  stress  (in  analogy  to  the  voltage)  to  the  displacement  velocity  (in  analogy  to  the 
current): 


Z(x1;  0)  = 


Tfm(x„0) 

V3m(xl5  0) 


(42) 


where  v£m  is  the  particle  displacement  velocity  in  the  2-2  composite.  Eq.  (42)  is  used  widely  in  the 
equivalent  circuit  model  of  transducers.36'38  Apparently,  for  a  2-2  composite  considered,  Z  from  eq. 
(42)  is  a  function  of  x1  due  to  the  variation  of  T3  and  v3  in  that  direction.  To  eliminate  this  variation,  the 
approach  taken  by  Miller  and  Pursey  is  adopted  here.33  That  is,  the  averaged  T3  and  averaged  v3  in  the 
x,-direction  are  used: 

7  _  JT3Cm(Xi,0)dxi 

M  M*(xlf  0)dx1  (43) 


where  the  integration  is  taken  in  one  unit  cell.  This  expression  seems  physically  meaningful  since  it 
reflects  an  averaged  mechanical  impedance  of  the  composite  at  the  interface.  Zin  calculated  from  eq. 
(43)  for  composites  with  44%  and  15%  ceramic  contents  is  shown  in  figure  20  which  is  quite  close  to 
those  in  figure  19.  In  the  following  discussion,  eq.  (43)  is  used  to  calculate  Zin  in  the  composite. 


One  interesting  question  one  might  ask  is  whether  the  acoustic  impedance  of  the  fluid  phase 
will  affect  the  effective  input  impedance  of  the  composite  at  the  interface.  To  answer  this  question,  the 
effective  input  impedance  of  the  composite  at  the  interface  is  evaluated  for  fluid  phases  with  different 
acoustic  impedance  (by  varying  the  bulk  modulus)  and  the  results  are  presented  in  figure  2 1 .  Clearly, 

the  effective  acoustic  impedance  of  a  2-2  composite  is  independent  of  the  acoustic  impedance  of  the 
fluid  at  the  interface. 


In  order  to  shed  light  on  the  large  change  of  the  input  acoustic  impedance  with  frequency,  the 
surface  displacement  u3  is  evaluated  at  the  center  of  the  ceramic  plate  (x,=0)  and  polymer  plate  (x,=d/2) 
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Figure  19.  (a)  The  effective  input  acoustic  impedance  Zin  of  the  2-2  composite  with  44%  ceramic 
content  calculated  from  the  data  in  figure  3  (the  reflection  coefficient).  For  the  comparison,  Zin 
calculated  from  eq.  (14)  where  only  the  amplitude  of  R  is  used.  The  results  demonstrate  the 

importance  of  the  phase  information  which  (b)  Zin  for  the  2-2  composite  with  15%  ceramic  volume 
content. 
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Figure  20.  The  effective  input  impedance  Zin  calculated  from  eq.  (16)  for  (a)  the  composite  of  44% 
ceramic  volume  content  and  (b)  the  composite  of  15%  ceramic  content.  The  results  here  are  nearly 
the  same  as  those  in  figure  5  except  at  frequencies  above  the  first  lateral  mode. 


Figure  21.  .  Zin  of  the  composite  of  44%  volume  content  derived  from  the  reflection  coefficient 
for  fluid  media  with  different  acoustic  impedance:  Zt  =  1.48  Mrayls  (water),  and  Z,=  2.09 
Mrayls.  Apparendy,  Zm  of  the  composite  for  the  two  cases  overlaps  with  each  other  indicating 
that  it  does  not  change  with  the  acoustic  impedance  of  the  fluid  medium. 


and  the  results  are  presented  in  figure  22(a).  At  low  frequencies,  the  surface  displacements  in  the 
polymer  and  ceramic  regions  are  in  unison  which  indicate  that  the  isostrain  model  used  in  many  earlier 
modelings  on  ultrasonic  composite  transducers  is  valid  in  this  frequency  region.3  30  As  the  frequency 
increases,  even  at  frequencies  still  far  below  the  first  lateral  mode  (at  fd  near  1),  the  surface  vibration 
amplitudes  in  the  two  regions  become  quite  different.  At  the  frequency  near  the  first  lateral  mode,  the 
vibrations  in  the  polymer  and  ceramic  regions  are  180  out  of  phase  and  in  this  frequency  region,  the 
effective  input  acoustic  impedance  of  the  composite  is  very  near  that  of  water  and  the  reflection 
coefficient  exhibits  a  minimum  as  shown  in  figure  17. 

As  we  have  pointed  out  in  the  introduction,  for  an  ultrasonic  piezoceramic  polymer  composite, 
one  of  the  most  important  factors  in  determining  the  performance  is  the  effective  acoustic  energy 
exchange  between  the  ceramic  and  the  external  medium.  As  we  have  seen  in  figure  17  and  19  where  at 
frequencies  near  the  lateral  mode,  the  input  acoustic  impedance  of  the  composite  is  close  to  that  of 
water  and  the  reflection  coefficient  reaches  a  minimum,  the  question  is  how  much  of  that  energy  enters 
into  the  ceramic  plates.  Figure  22(b)  shows  how  the  total  acoustic  powers  entering  into  the  ceramic 
region  and  polymer  region  vary  with  frequency.  Clearly,  at  the  high  frequency  region,  in  spite  of  the 
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fact  that  the  effective  transmission  coefficient  of  the  composite  increases,  the  amount  entering  into  the 
ceramic  plate  actually  decreases. 

It  is  also  interesting  to  examine  how  this  acoustic  energy  distribution  in  the  polymer  and 
ceramic  regions  changes  along  the  x3-direction  since  the  interaction  between  the  ceramic  and  polymer 
through  the  joint  region  will  cause  the  acoustic  energy  transfer  between  the  two.  In  figure  23,  the 
acoustic  energy  distribution  along  the  x3-direction  at  three  frequencies  (  f  «  f(,  f  ~  ft  /2,  and  f  —  f 
where  f,  is  the  first  lateral  mode  frequency)  is  shown.  In  all  the  cases,  there  is  a  redistribution  of  the 
acoustic  energy  along  the  x3-axis  and  the  acoustic  energy  in  the  polymer  is  gradually  transferred  to  the 
ceramic  plate.  However,  at  low  frequencies,  the  transition  region  is  much  shorter  than  that  at  high 
frequencies  (with  respect  to  the  wave  length  at  that  frequency).  For  instance,  in  figure  23(a),  the 
reduced  length  x3/X  is  about  0.01  for  the  acoustic  power  in  the  ceramic  to  reach  90%  of  the  final  value, 
while  in  figure  23(b),  this  region  increases  to  about  0.1,  and  in  figure  23(c),  it  becomes  near  0.2.  It 
should  be  reminded  that  in  a  thickness  mode  transducer,  the  thickness  of  the  piezocomposite  is  X/2. 

In  figure  24,  the  change  of  the  displacement  u3  in  the  polymer  and  ceramic  regions  along  the 
^■direction  at  the  three  frequencies  is  presented  and  it  shows  that  at  high  frequencies,  the  displacement 
amplitude  in  the  polymer  and  ceramic  regions  is  no  longer  the  same  even  deep  inside  the  composite. 

The  change  of  the  input  acoustic  impedance  from  eq.  (43)  along  the  x3-direction  in  the 
composite  is  also  evaluated.  Figure  25  are  the  results  of  the  2-2  composite  with  44%  ceramic  content  at 
different  frequencies.  Apparently,  at  high  frequencies  there  is  a  large  change  of  the  effective  acoustic 
impedance  from  the  surface  into  the  interior  of  the  composites  and  even  in  the  interior  of  the 

composite,  the  effective  impedance  is  not  the  same  as  that  derived  from  the  effective  medium 
theory.3'30 
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Figure  22.  (a)  The  surface  vibration  distribution  as  a  function  of  frequency  for  the  2-2  composite 
with  44%  ceramic  content  at  the  interface  with  water  and  under  an  acoustic  field.  u3p  and  u3c  are  the 
displacement  at  the  center  of  the  polymer  plate  (x,=d/2)  and  ceramic  plate  (x,=0).  (b)  The  acoustic 
energy  (in  reduced  units)  entering  the  ceramic  and  polymer  regions  as  a  function  of  frequency 
evaluated  at  the  interface  (x3=0)  where  P0  the  total  energy  in  the  incident  wave. 
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Figure  23.  The  energy  redistribution  inside  the  composite  along  the  x3-direction  where  x3=0  is  the 
interface.  There  is  an  energy  transfer  from  the  polymer  to  the  ceramic  inside  the  composite.  As 

shown  in  the  figures,  the  width  of  the  transition  region  increases  as  the  frequency  increases  (in 
reduced  length  unit). 
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Figure  24.  The  vibration  profile  as  a  function  of  the  distance  from  the  interface  with  water  (x3=0) 
where  u3p  and  u3c  are  the  displacement  at  the  center  of  the  polymer  plate  (x,=d/2)  and  ceramic  plate 
(xi=0).  At  low  frequencies,  the  polymer  and  ceramic  vibrate  with  the  same  amplitude  and  phase 
while  at  high  frequencies,  even  far  away  from  the  interface,  the  vibration  amplitude  of  the  two  is  still 
not  the  same  while  the  phase  becomes  the  same. 
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Figure  25.  The  effective  input  impedance  from  eq.  (16)  as  a  function  of  the  distance  from  the 
interface  with  water  (x3=0)  at  three  typical  frequencies:  far  below  the  lateral  mode,  at  about  the  half 
of  the  lateral  mode  frequency,  and  near  the  lateral  mode  frequency. 
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3.2.  Wave  reflection  and  transmission  at  solid-2-2  composite  interfaces: 


In  general,  a  piezoelectric  ultrasonic  transducer  has  a  multilayer  structure  in  which  there  are  one 
or  two  quarter  wavelength  impedance  matching  layers  in  front  of  the  piezoelectric  element  and  a 
backing  material  in  the  back.  The  piezoelectric  element  in  the  transducer  does  not  contact  directly  with 
fluid  medium.  Therefore,  the  investigation  of  the  interaction  between  a  2-2  composite  and  solid 
medium  has  practical  importance  for  composite  transducers. 

The  issues  investigated  in  this  section  are  quite  similar  to  those  studied  in  the  preceding  section. 
The  difference  is  that  a  solid  medium  can  support  shear  waves  while  a  fluid  cannot.  Due  to  this 

difference,  the  boundary  conditions  (3)  must  be  modified  to: 

Tl Cm  _  yS  yCm  _  yS  Cm  _  S  „Cm  S 
i  i33  ~  •‘33*  ■‘13  ~T13>U3  =«3’“l  =  «i 

4>c“  =  <(,*  £p  =  of ,  (44) 

where  the  superscript  S  refers  to  the  variables  in  the  solid  medium.  Hence,  the  variational  formula  for 
the  problem  becomes 

-  r3i“)(&3  +  Su£mf  +(7jS3  -  7j P)(&,p‘  +  &f  )* 

(45) 

+ (ufm  -  uixsr£m + st3s3)’ + (ufm  -  uf)(5r£m  + 

+(D 3S  -  Z>3Cm)(^5  +  S$Cm)*  +(<I>5  -  <Z>Cm)(6D$  +  8D%mf}dS.  =  0 

where  the  surface  integral  is  performed  at  x3=0  over  one  repeat  unit  in  the  x,-direction  of  the 

composite.  The  elastic  displacement  wf,wf,  and  the  electrical  potential  d^in  the  solid  medium  are 

expanded  in  terms  of  the  eigenfuntions: 

2 

ul  =  jl  1 8ni  sin {h%x{)  exp(-y^.r3  )R% 

n  i= 1 

2  (46> 
u3  =  Po  expO'/5ox3)-;I  I  /J  cos(h%Xl)  cxp(-jp^x3)R^ 

n  z=l 

d>5  =  lcos(hSnx{)exp(h%x3)R% 
n 

where  exp(jpQX3)  is  the  incident  wave  and  other  terms  are  from  the  reflected  waves  in  which  both 
the  longitudinal  and  shear  components  exist  except  the  one  with  n=0  where  the  reflected  wave  is  a  pure 
longitudinal  one.  Substituting  the  elastic  and  electrical  variables  in  the  2-2  composite  and  eq.  (46)  into 
equation  (45),  the  coefficients  of  the  expansions  can  be  obtained. 

Figure  26  presents  the  dispersion  curve  of  the  reflection  coefficient  at  the  interface  for  2-2 
piezocomposites  with  44%  ceramic  volume  content  and  the  acoustic  impedance  calculated  from  eq. 
(41)  where  the  characteristic  impedance  of  the  solid  medium  is  Z=2.4  Mrayls. 
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By  comparing  the  results  in  figure  26  and  figure  19,  it  can  be  found  that  the  input  acoustic 
impedance  of  2-2  composites  calculated  here  is  different  from  that  in  the  fluid  case.  In  addition,  it  is 
also  found  that  at  a  solid  medium-2-2  composite  interface,  the  input  acoustic  impedance  also  depends 
on  the  characteristic  impedance  of  the  solid  medium  at  the  interface,  which  seems  to  be  quite  different 
from  the  fluid  medium  interface. 

To  understand  this  phenomenon,  we  notice  that  the  difference  between  a  solid  medium  and  a 
fluid  one  is  that  a  solid  medium  has  a  non-zero  shear  stiffness  constant.  Therefore,  Zin  for  a  2-2 
composite  (44%  ceramic  content)  in  contact  with  solid  medium  is  evaluated  where  p  (mass  density) 
and  cu  are  kept  as  constant  and  c^  is  varied  and  the  results  are  presented  in  figure  27(a).  Similarly 
Zin  f°r  a  2-2  composite  in  contact  with  solid  medium  is  also  evaluated  where  p  and  c44  are  kept 
constant  and  cu  is  varied  and  the  results  are  shown  in  figure  27(b).  The  results  demonstrate  that  for  a 
plane  incident  wave,  even  if  the  characteristic  longitudinal  impedance  of  the  solid  medium  is  kept 
constant,  the  effective  input  impedance  of  a  composite  and,  hence,  the  reflection  coefficient  from  the 
interface  will  change  if  the  shear  stiffness  coefficient  of  the  solid  medium  changes.  On  the  other  hand, 
if  the  shear  stiffness  constant  is  kept  constant  and  the  characteristic  longitudinal  impedance  is  varied  in 
the  solid  medium,  the  effective  input  acoustic  impedance  of  a  composite  will  not  have  much  change. 
This  phenomenon  is  directly  related  to  the  surface  uniformity  of  the  composite  and  as  has  been  shown 

earlier,  the  enhanced  stress  transfer  between  the  ceramic  and  polymer  regions  is  through  the  shear 
action.39 


To  illustrate  this,  the  surface  vibration  distribution  in  the  ceramic  and  polymer  regions  is  also 
evaluated  for  the  2-2  composite  in  contact  with  solid  media  of  different  cu  and  c^  and  the  results  are 
presented  in  figure  28.  There  is  veiy  little  difference  in  the  surface  vibration  profile  for  solid  media 

with  the  same  c^  and  different  c„.  In  contrast,  the  surface  vibration  profile  will  change  as  c^  is 
changed. 

Figure  29  illustrates  the  input  acoustic  impedance  of  a  composite  as  a  function  of  x3  for 
different  solid  media  at  the  interface.  It  shows  the  length  of  the  transition  region  in  the  composite 
where  the  evanescent  waves  are  important  decreases  as  the  acoustic  impedance  of  the  medium, 
especially  the  shear  modulus,  increases. 
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f*d  (MHz* mm) 


Figure  26.  (a)  The  reflection  coefficient  at  a  interface  between  a  solid  medium  (p=  ,  cu=  c  =  ) 

and  a  2-2  composite  of  44%  ceramic  volume  content  as  a  function  of  frequency,  (b)  The  effective 
input  acoustic  imepdance  of  a  2-2  composite  of  44%  ceramic  volume  content  as  a  function  of 
frequency.  Z  is  calculated  from  the  reflection  coefficient. 


Figure  27.  Data  illustrate  the  dependence  of  the  input  acoustic  impedance  of  a  composite  on  the 
elastic  properties  of  the  solid  medium  at  the  interface,  (a)  The  solid  medium  for  the  curve  1  is 
p-1.61  g/cm  ,  cu  =  1.36*  1010  N/m2,  and  C+4  =  3.43*  109  N/m2,  and  for  the  curve  2  and  curve  3 
both  p  and  cn  are  kept  constant  and  cw  =  6.86*109  N/m2  and  c44  =  1.72*109  N/m2.  (b)  Here  both 
p  and  C44  are  kept  constant,  and  ctl=  0.68  *1010,  cu=  1.36  *1010,  and  cn=  2.72  *10ION/m2.  Hie 
amplitude  of  Zjn  does  not  change  much  with  c,,  except  a  small  change  in  the  phase. 
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FIr.f1 28  V  Ihe  flguresmustrate  the  effect  of  cn  and  Cj4  of  the  solid  medium  on  the  surface  vibratiot 

P  .  a',°  the3comP“1te-  (a)  The  parameters  of  the  solid  medium  for  the  data  of  dashed  curve  is 

P_1.6  g/cm  •  cu  =  13610  N/m  ,  and  cu  =  3.43*10’  N/m!,  and  for  the  solid  curve,  both  p,  c 

“rC  *  S,™e.^„C?pt  cii-232*10  ■v/rn3  (b)  The  parameters  for  the  solid  medium  is  p=1.61  g/cmt 

om mJ rf,n10  N/m''  and  C"  =  6'86*10  N/mI-  Noti«  affect  of  Cj4  on  the  surface  vibradon 
prome  or  the  composite. 


Figure  29.  The  figures  illustrate  the  effect  of  the  elastic  properties  of  the  solid  medium  at  the 
interface  on  the  effective  input  impedance  distribution  along  the  x3-axis  at  a  frequency  f*d  =  0.446 
Mhz*mm  where  x3=0  is  the  interface,  (a)  The  dashed  line:  p=1.16  g/cm3,  cu=  7.72*109  N/m2,  and 
044  =1.588*109  N/m2.  The  solid  line:  p=l.l  g/cm3,  cn  =  5.41*109  N/m2,  and  c^  =  1.307*109  N/m2 
(Spurrs  epoxy),  (b)  The  parameters  for  the  solid  medium  is:  p=1.61  g/cm3,  ctl  =  1  36*1010  N/m2 
and  C44  =  3.43*109  N/m2. 


3A  ..Reflection  and  transmission  from  the  boundary  between  fluid  and  acoustic  impedanrp 
matching  layer: 


We  now  proceed  to  investigate  how  the  various  effects  observed  in  the  preceding  sections 
affect  the  selection  of  anti-reflection  matching  layers  at  the  interface.  Figure  30  is  a  schematic  of  an 
acoustic  system  which  consists  of  a  fluid  medium,  a  matching  layer  and  a  2-2  composite.  In  this  case, 
there  are  two  interfaces  with  different  boundary  conditions  to  be  treated.  One  of  the  interfaces  is 
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Figure  30.  Schematic  of  an  interface  system  consisting  of  a  fluid  medium  (medium  1) 
matching  layer,  and  a  2-2  composite.  The  thickness  of  the  matching  layer  is  t  and  the 
period  of  the  composite  is  d. 

between  the  2-2  composite  and  matching  layer  where  the  boundary  conditions  are  those  of  equation 
(44).  Another  one  is  the  fluid  and  matching  layer  interface  where  the  boundary  conditions  are: 

T33  =  33 . 713  =  0,  «3S  =  =  D?  at  x3=-t  (47) 

where  t  is  the  thickness  of  the  matching  layer  and  the  other  notations  are  the  same  as  before.  The 
variational  formula  for  this  problem  can  be  derived: 

is,  f  + &f")* + (Tii  -  7if'”x&p‘ + &-y 

+(“3°"  -4  xsiS" +«$>* + (“fm  -  “h<srxSm + sr&f 

+ (Df  -  a ?”x»s + »°”)*  +(os  -  + 8$")'ys  (48) 

+is,  21<TH  ~  r33X&3  +  &?')’  +  IT^iSuf +  (uf  -  )(Slf3  +  sr£ )* 

+  (D?  -  D$ )(<5G>5  +  <5$^)*  +  (®w  -  ®s)(SD$  +  5D^)*]dS  =  0 
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where  the  surface  integrals  Sl  and  S2  are  performed  at  x3=0  and  x3=-t  planes,  respectively,  from  Xj=0 
to  Xj=d/2. 


In  the  matching  layer  region,  there  are  waves  which  propagate  in  both  the  +x3  and  -x3 
directions  due  to  the  two  interfaces  at  x3=0  and  x3=t.  Hence,  the  expressions  for  the  elastic  and 

electrical  variables  become: 

2 

U1  =JZ  I gni (Pni ) sin^n  *1 )[- txp(-j/3^x3 ) Rln  +exp(y^x3)^] 
n  i=l 

2 

«3  =  I  X  fniipni ) cos(^f jq ) [- exp(-y^x3 ) RXn  +  txpij^x^R^]  (49) 

n  i~l 

Os  =  Icos(/^x1)[-exp(-/^x3)^+exp(^fx3)i?2] 

n 

*The  solutions  in  the  composite  and  in  the  fluid  medium  are  those  of  eqs.  (37),  (38)  and  (39). 
Substituting  those  expressions  into  equation  (48),  the  undetermined  coefficients  in  these  equations  can 
be  found  and  various  quantities  in  different  regions  can  be  evaluated. 


In  general,  for  homogeneous  materials,  the  reflection  coefficient  from  the  boundary  between  a 
fluid  and  a  matching  layer  will  be  zero  when  the  thickness  of  the  matching  layer  is  quarter  wavelength 
and  its  acoustic  impedance  ZM  is 


zM  =  Vz^F 


(50) 


where  Z  and  Z  are  the  acoustic  impedance  of  the  media  on  the  two  sides  of  the  matching  layer, 
respectively.19,31 


For  a  composite  material,  from  the  results  obtained  in  the  preceding  sections,  it  is  expected  that 
the  parameters  of  the  anti-reflection  matching  layer  will  be  modified.  Figure  31  illustrates  the  derived 
reflection  coefficient  from  the  boundary  between  the  fluid  and  matching  layer  for  different  matching 
layer  materials  where  the  thickness  of  the  matching  layer  is  chose  so  that  the  frequency  of  the  reflection 
coefficient  minimum  is  at  fd=0.335  MHz*mm.  From  results,  it  is  found  that  the  acoustic  impedance  of 
the  matching  layer  from  which  the  reflection  coefficient  becomes  zero  satisfies  approximately  equation 
(50)  if  Zs  is  replaced  by  the  effective  input  acoustic  impedance  of  the  2-2  composite  at  the  low 
frequency.  However,  the  thickness  of  the  matching  layer  is  no  longer  quarter  wavelength  and  in  fact,  it 
is  shorter  than  774.  These  results  are  summarized  in  figure  32.  The  data  in  figure  32(a)  is  the  acoustic 
impedance  of  the  matching  layer  in  order  for  the  reflection  coefficient  to  be  zero  as  a  function  of 
frequency.  In  figure  32(b),  the  change  of  the  matching  layer  thickness  as  a  function  of  frequency  is 
presented  where  t=l  corresponds  to  the  thickness  of  774  at  that  frequency.  Therefore,  as  the  frequency 
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Pd  (MHz* mm) 


Figure  31.  The  reflection  coefficient  from  the  matching  layer  system  where  the  total  transmission 
occurs  at  f*d  =  0.334  MHz*mm  for  the  matching  layer  with  t  =  2.02*d  =  0.91*(A74)  (at  f*d  =  0.334 
MHz*mm)  and  ZM=4.68  Mrayls  (p=1.61  g/ cm3,  cu  =  1.36*1010  N/m2,  and  cw  =  3.43*109  N/m2) 
(the  curve  1).  For  the  comparison,  the  reflection  coefficient  from  the  matching  layers  with  other 
parameters  is  also  shown  here:  for  the  curve  2,  ZM  =  3.0  Mrayls  and  t  =  1.85*d,  and  for  the  curve 
3,  ZM  =  2.44  Mrayls  and  t  =  1.58*d. 

increases,  both  the  acoustic  impedance  and  the  thickness  of  the  matching  layer  decrease.  For  example, 
when  the  frequency  changes  from  0.05  MHz*mm  to  0.5  MHz*mm,  the  acoustic  impedance  of  the 
matching  layer  changes  from  4.70  Mrayls  to  4.42  Mrayls  and  the  thickness  t  changes  from  about  A/4 
to  0.86*A/4.  The  reduction  in  the  matching  layer  thickness  is  a  direct  result  of  the  fact  that  the  incident 
plane  wave  suffers  more  than  180°  phase  loss  upon  the  reflection  from  the  interface. 

From  the  fact  that  in  the  frequency  range  of  interesting,  the  evanescent  waves  at  the  interface 
wifi  decay  approximately  as  exp(-2mcx3/d),  the  effect  of  the  interface  of  x3=t,  where  t  is  approximately 
A/4,  to  the  surface  vibration  distribution  at  the  composite-matching  layer  interface  will  be  quite  small. 
Therefore,  the  reflection  coefficient  introduced  in  the  section  IE  between  the  composite  and  a  solid 
medium  can  be  used  to  approximate  the  reflection  at  the  matching  layer-composite  interface.  Hence, 
the  reflection  coefficient  R  from  the  system  in  figure  30  can  be  approximated  as: 

D_  R^  +  R^e2^ 

l+R^e2^*  -  (51) 

where  R12  and  R23  are  the  reflection  coefficients  at  the  fluid-matching  layer  interface  and  matching 
layer-2-2  composite  interface,  respectively.  In  eq.  (5 1 ),  R12  is  the  reflection  coefficient  from  the 
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Figure  32.  The  required  parameters  for  (a)  ZM  and  (b)  t  of  the  matching  layer  to  have  total 
transmission  as  a  function  of  frequency  for  a  water-2-2  composite  (44%  ceramic  volume  content) 
interface. 


interface  of  two  semi-infinite  media  of  the  fluid  and  solid,  and  R23  is  that  derived  in  the  section  III  (the 
reflection  coefficient  from  the  solid  medium-composite  interface).  From  equation  (51),  the 
requirements  to  the  matching  layer  to  achieve  the  total  transmission  (R=0)  can  be  derived.  The  acoustic 
impedance  of  the  matching  layer  is 


zwlzCml 

zCm 

-Zw  cos(0) 

zCm 

cos(6)-Zw 

(52) 


where  6  is  the  phase  angle  of  the  acoustic  impedance  of  the  2-2  composite.  The  thickness  of  the 
matching  layer  is 


Xn  1 

— {1  — arctan( 
4  K 


2Zm 

zCm 

sin(i3) 

zCm 

2-(zM)2 ' 

■)} 


(53) 


iAs  shown  in  figure  32,  in  the  frequency  range  of  interest,  the  results  from  eqs.  (52)  and  (53)  are 
almost  identical  to  those  derived  from  the  numerical  calculation. 


TV.  Modeling  and  Experimental  Study  of  1-3  Piezocomgosjtes: 


Both  analytical  modeling  and  experimental  investigation  of  1-3  composites  are  currently 
underway.  As  a  first  step  to  model  a  1-3  composite,  a  concentric  unit  cell  as  schematically  shown  in 
figure  33  is  adopted.  This  unit  cell  is  an  approximation  to  a  1-3  composite  with  the  ceramic  rods 


arranged  in  a  hexagonal  lattice.  In  the  analysis  here,  the  elastic  loss  in  the  polymer  phase  is  included. 
The  other  losses  such  as  the  dielectric  loss  in  the  ceramic  phase  can  also  be  included.  However. 


compared  with  the  elastic  loss  in  the  polymer  phase,  these  other  loss  mechanisms  may  not  be 
significant.  The  wave  equation  and  the  boundary  conditions  can  be  expressed  as 

pur  =  Trrr  +  TrZtZ  +  - 


and 


Trr  ~  Tqq 


P^z  ~  Trz,r  + 

Dnr  +  ?L 


irz+r 

r  +1zz.z 


(54) 


+  D7  r  —  0 


jC  —  ,tP  tC  _  -t -P  ff.C  (f\P 
irr  -  irr,  irz  -  ir.,q)  =  <p  , 

C  P  C  P  C  P 
ur  =uf,u^  =  uC,  D}z  =Df, 

D, r  =  0,  uf  =  0,  TrPz  =  0, 


at  r=r, 


at  r=r2 


(55) 
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The  constitutive  relations  for  ceramic  phase  are 

TrCr  =  cfi«£r  +  4  £  +  -  e3l  E? 

T0d  =  cnur,r  +  cn  L-y  +  c13uZz  ~  *3 1  Ez 

c 

TS  =  cl3^ur,r  +  ~y)  +  c33uZz  ~  e33 Ez 

=  c44(ur,z  +  uz,r ) "  el5Er 

Dr  =  *15  (“£z  +  «£■)  +  4\EZ 

c 

Dz  =  *31  (uZr  +  ~)  +  *33  uz,z  +  £33ez 


(56) 


The  constitutive  relations  of  the  polymer  phase  are  obtained  by  setting  the  piezoelectric  coeficients 

e3i’e33>  30(1  eis  zero  in  eq.  (39).  The  superscripts  C  and  P  denote  the  ceramic  and  polymer  phases, 
respectively. 


From  the  symmetry  consideration,  the  solutions  to  the  equation  (54)  are 
u9  =  A7 j  (hEr)  cos(/fe)  exp(jcot) 

uz  =  BJQ(h^r)sin(Pz)  exp(jcot)  (57) 

=  CJq  (hCr)  sin(Pz)  exp(jcot) 

where  A,  B,  and  C  are  constants,  J0  and  J,  are  zeroth  and  the  first  order  Bessel  functions  respectively, 
h  and  (3  are  wave  vectors  along  the  radial  and  the  axial  directions,  respectively. 

Substitude  Eq.  (56)  and  (57)  into  Eq.  (54)  yields  three  homogenieous  equations  which,  for  a 
given  co  and  p,  yield  three  hc,  corresponding  to  the  quasi-longitudinal,  quasi-shear,  and  quasi- 

electromagnedc  waves  in  the  ceramic  phase.  Hence,  the  general  solutions  in  ceramic  phase  can  be 
expressed  as 

2 

ur  =  I  CtfiJ j  (hfr)  cos(fiz) 
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(58) 


3 

uz  -  £  CifiJo(hf  r)sin(/fe) 
i=l 

3 

=  X  CftjQihf r)sin(pz) 
i'=l 

For  the  polymer  phase,  the  solutions  are 
2 

ur  ~  X  J\(hf  r)  +  Bjgf  Y^(hf  r)]cos(j3z) 

i'=l 
2 

uz  =  llAif  Jo(hfr)  +  BifiPY0(hf, r)]s  in(ft)  (59) 

^  =[^o(^r)  +  ^Jo(^3>r)]sin(/fe) 

where  Y0and  Y,  are  modified  zeroth  and  the  first  order  Bessel  functions  and  hp  are 


h/=j(3 


The  elastic  stiffness  of  polymer  can  be  expressed  as 
cuP=  cu-jTinco 
Cu?=  C^-jTl^CO 

where  Tjjj,  are  the  viscosity  constants,  by  the  boundary  conditions  (eq.  (55)). 

Substituting  Eq.  (58)  and  Eq.  (59)  into  the  boundary  conditions  (55)  yields  9  homogenieous 
equations  from  which  the  relations  among  the  coefficients  Ai5  B;,  D;,  and  C,  can  be  determined  which 
are  functions  of©,  (3  and  the  material  parameters  of  both  ceramic  and  polymer  phases: 

(Ky)  (Aj)  =  0 

where  (Kij)  is  a  9x9  matrix  and  (Aj)  a  1x9  matrix  in  which  the  elements  are  C,,  C2,  C3,  A,,  A2,  B,, 
®2’  ^2.  The  condition  for  a  non-trivial  solution  to  these  equation  is 

IKyl  =  0, 

from  which,  the  relation  between  (3  and  ©,  i.e.,  the  dispersion  curves  for  the  guided  waves  in  a 
composite,  can  be  determined.  These  procedures  are  exactly  parallel  to  those  in  a  2-2  composite. 
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Shown  in  figure  34  are  the  dispersion  curves  for  a  1-3  composite  with  44%  ceramic  volume 
content  where  the  loss  in  the  polymer  phase  is  assumed  zero  (PZT-5H  and  Spur  epoxy).  However, 
the  dispersion  curves,  especially  the  ones  related  to  the  periodic  structure  in  the  lateral  direction,  show 
changes  with  the  elastic  loss.  As  presented  in  figures  35  and  36,  the  cut-off  frequency  for  the 
dispersion  cuver  for  the  second  branch  (the  lateral  mode)  becomes  zero  as  Pd  approaches  zero  when 
the  elastic  loss  in  introduced.  On  the  other  hand,  the  imaginary  part  of  the  wave  vector  increases 
rapidly  as  (3d  -»  0,  indicating  that  in  the  limit  (3d  ->  Othe  second  mode  (the  lateral  mode)  becomes 
non-propagational. 

Experimentally,  it  was  observed  that  the  behavior  of  the  lateral  modes  in  a  1-3  composite  is 
quite  different  from  the  prediction  based  on  the  stop-band  edge  resonance  and  seems  to  be  much  more 
complicated  than  that  in  a  2-2  composite.  These  results  will  be  presented  in  the  future. 


▲ 


Figure  33  Schematic  of  a  1-3  composite  which  approximates  the  unit  cell  of  a  1-3  composite  with  a 
hexagonal  lattice. 
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jure  36.  The  effect  of  increased  elastic  loss  (compared  with  figure  35)  in  the  polymer  phase  on  the 
dispersion  curves  where  the  parameters  for  the  polymer  phase  are:  p=1100  kg/m3,  c.,=  5.41xl09 
N/m  ,  C44=  1.037xl09  N/nr,  riu=  101.12  N  s/m2,  ri^=  36  N  s/m2 


*  >  r 


V.  Summary: 

A  dynamic  model  on  piezoceramic  polymer  composites  with  laminar  periodic  structure  are 
presented.  The  result  shows  that  the  various  resonant  modes  in  a  composite  stmcture  can  be  traced 
back  to  the  modes  in  either  an  isolated  ceramic  plate  or  polymer  plate  with  appropriate  boundary 
conditions  (stress  free  for  the  ceramic  plate  and  strain  free  for  the  polymer  plate).  It  also  shows  that 
there  exist  a  series  of  modes  associated  with  the  periodic  structure  of  a  composite,  which  is  beyond  the 
stop-band  edge  resonance  prediction.  One  of  the  main  concerns  in  designing  a  composite  transducer  is 
how  the  surface  vibration  profile  changes  with  frequency  and  how  this  is  influenced  by  the  aspect  ratio 
t/d.  It  was  predicted  and  verified  by  experiment  that  as  long  as  the  thickness  resonance  is  below  the 
first  lateral  mode  frequency,  there  is  always  a  frequency  f,  which  is  near  the  thickness  resonance  and  at 
..which  the  polymer  and  ceramic  vibrate  in  unison.  The  effect  oft/d  is  to  change  the  position  of  f,  with 
respect  to  the  thickness  resonance  frequency  and  the  bandwidth  in  which  polymer  and  ceramic  have 
nearly  the  same  vibration  amplitude  and  phase.  It  is  also  predicted  that  when  operated  in  a  fluid 
medium  such  as  water,  there  will  be  a  resonance  mode  whose  frequency  is  determined  by  the  velocity 
of  the  fluid  medium  and  the  unit  cell  length  d  and  is  associated  with  the  oscillation  of  the  fluid,  causing 
the  polymer  and  ceramic  to  vibrate  180°  out  of  phase.  The  difference  in  the  surface  vibration  profile! 
between  in  air  and  in  water  indicates  the  need  to  characterize  the  vibration  pattern  of  a  composite  in  a 
fluid  medium  since  it  is  much  closer  to  the  real  application  environment. 

The  reflection  and  transmission  of  a  plane  acoustic  wave  at  a  medium-composite  interface  and 
the  issues  related  to  the  design  of  matching  layer  for  a  composite  are  analyzed  based  on  the  model 
developed.  It  has  been  shown  that  the  reflection  coefficient  from  the  interface  is  a  complex  and  the 
reflected  wave  suffers  a  more  than  180°  phase  change.  The  effective  input  acoustic  impedance  Zin  of 
the  composite  at  the  interface  was  evaluated  and  both  the  amplitude  and  phase  show  a  strong  frequency 
dependence.  For  a  fluid  medium,  it  was  found  that  Zm  does  not  change  with  the  acoustic  impedance  of 
the  medium.  However,  for  a  solid  medium,  Zin  will  change  if  the  shear  stiffness  constant  of  the 
medium  changes.  It  was  demonstrated  that  this  difference  originates  from  the  non-uniformity  of  the 

surface  vibration  distribution  of  the  composite  at  the  interface  which  depends  crucially  on  the  shear 
stiffness  constant  of  the  medium. 

Since  for  a  piezoceramic  polymer  composite,  it  is  the  ceramic  phase  which  performs  the  energy 
conversion  between  the  acoustic  and  electric  forms,  how  much  acoustic  energy  can  enter  the  cerandc 
region  is  one  of  the  most  important  parameters  in  a  composite  transducer  design.  In  the  paper,  we 
show  that  even  though  the  effective  transmission  coefficient  increases  as  the  frequency  is  increased, 
the  amount  of  acoustic  energy  entering  the  ceramic  region  actually  decreases.  Therefore,  there  may  be 
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a  trade-off  between  the  bandwidth,  which  is  related  to  the  transmission  coefficient,  and  the  sensitivity 
in  the  composite  transducer  design. 

From  the  fact  that  there  is  more  than  180°  phase  change  in  the  reflection  from  the  medium- 
composite  interface,  it  is  shown  that  the  matching  layer  thickness  is  no  longer  equal  to  the  quarter 
wavelength  but  smaller  than  that.  In  addition,  the  acoustic  impedance  of  the  matching  layer  will  also 
be  affected  by  the  phase  of  Zin  of  the  composite. 

1-3  piezocomposite  was  also  modeled  in  which  the  elastic  loss  of  the  polymer  phase  is 
included.  It  was  observed  that  at  small  Pd  limit,  the  cut-off  frequency  of  the  second  branch  which  is 
related  to  the  lateral  mode  approaches  zero  and  the  imaginary  part  of  the  wave  vector  increases  rapidly. 
^Experimentally,  it  was  observed  that  the  modes  related  to  the  periodic  structure  of  a  1-3  composite 
seems  to  be  much  more  complicated  than  the  prediction  of  the  stop-band  edge  resonance. 
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